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tions from  Mr.  R.  Smith  of  Shaker  Research.  Failure  data  collection  was 
assisted  by  Mr.  A.  Fries  of  the  Rotron  Manufacturing  Company.  The  RADC 
Project  Engineer  was  Mr.  L.  Gubbins. 
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EVALUATION 


This  contractual  effort  is  part  of  the  broad  RADC  Reliability 
Program  intended  to  provide  reliability  prediction  procedures  for 
military  electronic  equipment  and  systems.  These  prediction  procedures 
are  contained  in  MIL-HD6K-21 7B  for  which  RADC  is  the  preparing  activity. 
The  final  model  developed  in  this  study  will  replace  the  present 
electric  motor  reliability  models  that  are  now  in  MIL-H06K-217B. 

a U V ^ ^1**^ 

LESTER  J.  GUBBINS 
Project  Engineer 


1.0  INTRODUCTION 


llils  report  sun«narir,en  the  renultn  of  a study  directed  towards  the  de* 
velopment  of  a failure  rate  prediction  model  for  fractional  horsepower  electric 
motors.  Ihe  objective  of  the  study  was  to  completely  revise  the  Information 
and  failure  rate  models  for  fractional  horsepower  electric  motors  an  used  In 
electronic  equipment  currently  Included  In  Section  2.11.1  of  MIL-HDBK«217B 
dated  20  September  1974. 

To  accomplish  the  program  objectives,  an  existing  failure  data  bank 
(Reference  1)  of  over  2(X)0  failures  of  fractional  horsepower  motors,  was  up- 
dated and  supplemented  from  a survey  of  electric  motor  manufacturers.  The 
major  portion  of  the  failure  data  history  was  furnished  by  the  Rotron  Manu- 
facturing Company.  A regression  analysis  was  conducted  on  all  failures 
utilizing  such  variables  as  motor  power,  speed,  temperature,  failure  mode, 
Insulation  type  and  motor  type  to  establish  the  dominant  motor  life  Influencing 
parameters  and  their  Interrelationships. 

Two  competing  failure  modes  were  found  to  dominate  motor  life;  l.e,, 
bearing  failures  and  winding  failures.  All  bearings  In  the  failure  bank  were 
of  the  greace-packed  ball  bearing  type.  Separation  of  failure  by  failure  mode 
Indicated  that  both  fafure  modes  could  be  described  by  a Welbull  Commulative 
Distribution.  A mathematical  model  for  each  failure  mode  wl.th  dominant  varia- 
bles was  developed  and  a competing  rli:k  failure  model  established. 

A literature  aur*^  of  bearing  failures  and  winding  failures  was  con- 
ducted to  verify  failure  nodes  and  dependent  failure  variables.  The  follow- 
ing sections  of  this  report  describe  the  data  failure  bank,  failure  modes, 
literature  review  and  developnant  of  the  mathematlca'.  prediction  models.  Ap- 
pendix A describes  the  regression  analysis,  end  Appendix  B outlines  the 
recosBMnded  replseement  for  Section  2.8.1  of  MIL-HDBK-217B. 
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2.0  SUMMARY 


At,  existing  data  failure  bank  i-onsisting  of  over  2000  actual  fractional 
horsepower  (FHP)  motor  failures  was  upgraded  and  utilired  to  develop  a failure 
model  for  FHP  motors.  Motor  failures  were  dominated  by  two  failure  inodes; 
i.w. , winding  failures  and  bearing  failures. 

A mathematical  reliability  model  was  developed  for  each  failure  mode 
utilizing  a Weibull  cumulative  distribution  function  and  regression  technique. 
Results  of  a literature  survey  of  both  failure  modes  were  used  to  supplement 
failure  data  where  required  and  verify  the  developed  models. 

Failure  modes  %rare  recombined  into  a competing  risk  cumulative  distri- 
bution model  for  use  in  predicting  overall  motor  life  and  failure  rates. 
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3,0  DISCUSSION 


The  development  of  a "Reltabllity  Model"  for  fractional  horsepower 
(FHP)  ntctcrs  necessitates  « well “documented  data  bank  of  motor  failures  on 
which  to  base  failure  predictions.  To  develop  such  a data  bank,  47  manufac- 
turers of  fractional  horsepower  motors  were  contacted  for  assistance.  In 
addition,  a literature  search  was  conducted  through  the  Engineering  Index 
(El)  and  the  National  Technical  Index  Service  (NTIS)  covering  motors,  wind- 
ings, insulation,  bearings  and  failure  modes.  As  a result  of  these  efforts 
over  3000  failure  histories  were  collected  with  1614  failures  suitably 
documented  for  analysis  purposes.  The  major  portion  of  the  failure  bank 
was  furnished  by  the  Rotron  Manufacturing  Company. 

The  composition  of  motor  types  encompassed  by  the  failure  bank  covers 
power  ranges  from  15  watts  input  to  625  watts  input.  Motor  types  encompassed: 

Synchronous  polyphase  motors  Shaded  pole  motors 

Polyphase  induction  motors  Split  phase  motors 

Capacitor-run  motors  Capacitor-start  motors 

All  motors  analyzed  contained  grease  packed,  single  and  double  shielded 
ball  bearings  to  Antifriction  Bearing  Engineers  Committee  (ABEC)  quality  3 
and  5.  Motor  (NEMA)  insulation  classes  F and  H corresponding  to  IEEE  classes 
150°C  and  180^C  are  the  only  classes  used  in  the  failure  bank. 

Operating  parameters  of  the  motors  cover  the  steady  state  operating 
range  from  18**C  .n  140°C  and  thermal  cycling  from  -65°C  to  125**C,  Operating 
speeds  varied  Viom  1750  RPM  to  23000  RPM  covering  two  input  frequencies;  l.e., 
60  Hz  and  400  Hz.  Motor  input  voltage  Included  115V,  208V,  230V,  416V  and 
440V. 

Because  of  the  larc'  nuid>er  of  variables,  all  data  was  stored  in  a 
digital  computer  with  the  following  information  available  for  each  test 
point: 
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• Failure  mode 

• Time  to  failure 

• Date  of  failure 

• Motor  type 

• Bearing  class 

• Searing  lubricant 

• Insulation  class 


• Ambient  temperature 

• Bearing  temper'^ture  rlsa 

• Motor  temperature  rise 

• Motor  load 

• Motor  frame  size 

• Special  construction  feature 


3.1  Data  Bank  Analysis 

The  major  portion  of  the  failure  histories  had  previously  been  analyzed 
(Reference  1)  for  failure  modes  with  the  following  results. 

Bearing  failures  80.85% 

Electrical  failures  16.55% 

Mechanical  failures  2.60% 


These  results  clearly  Indicated  that  a failure  model  for  electric  motors 
could  be  developed  utilizing  two  basic  failure  modes;  l.e.,  bearing  failures  < 

and  winding  failutwi.  In  a normal  population  of  motors,  it  would  be  antici- 
pated that  the  dominant  failure  mode  would  be  the  result  of  bearing  failures. 

This  condition,  however,  Is  limited  to  non-regreaseable  ball  bearing  motors. 

As  a result  of  these  considerations  only  bearing  and  electrical  failures  were 
analyzed  from  the  data  bank. 

Failure  data  are  dominantly  composed  of  life  test  results  of  individual 

j 

teat  populations.  It  was  found  (1)  that  the  most  effective  way  of  handling 
the  data  was  by  means  of  a Welbull  Cumulative  Distribution  analysis  of  each 
individual  test  population.  The  results  of  this  analysis  provides  a linear 

regression  best  fit  Welbull  slope  (8)  and  characteristic  life  (a)  for  each  j 

% 

test  group  of  motors.  Additional  regression  techniques  are  applied  to  deter-  j 

mine  the  ''nfluence  of  parameters  such  as  temperature,  speed,  bearing  lubricant,  i 

motor  type,  etc.  on  characteristic  life.  I 


tion  from  test  groups  Influences  ^he  order  ranking  of  failures  inducing  errors 
in  the  predictions  of  Weibull  slope  and  characteristic  life.  A technique  for 
separating  failure  modes  as  proposed  by  Nelson  (2)  is  tc  consider  all  failures 
of  one  failure  mode  as  suspended  (censored)  test  units  vhen  analyzing  a 
second  failure  mode.  This  is  accomplished  by  means  of  a Weibull  Cunulatlve 
Hazard  function.  Censoring  of  failures  of  all  other  failure  modes  except  the 
failure  mode  of  interest  was  applied  to  each  test  group  population  during  the 
data  analysis.  Tlie  reader  is  referred  to  Reference  2 for  a comprehensive  dis- 
cussion of  censoring  and  only  the  definition  of  the  cumulative  hazard  function 
is  provided;  l.e. 

Weibull  Cumulative  Distribution  Function  is  defined  as: 


F(t)  “ 1 - e 


-(feP 


Probability  Density  Function  is  defined  as: 


f(t)  - ^ (F(t))  - ^ t ^"^e 


Hazard  Function  (failure  rate)  is  defined  as: 


l-F(t)  5ip  ^ 


Cumulative  Hazard  Function  is  defined  as: 


H(t)  - r t»(t)dt  - (^P 

Jri 


A further  data  separation  was  found  desirable  lu  tabulation  of  the  re- 
sults of  the  analysis;  l.e.,  tabulation  of  bearing  failures  by  grease  type. 
Three  basic  greases  were  utilised  for  the  majority  of  the  tests.  Approxi- 
mately 30  addltonal  grease  types  were  tested  In  motors,  but  the  majority  of  the 
greases  did  not  have  sufficient  test  points  to  provide  statistically  signifi- 
cant results. 
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Tables  i through  4 sunmarize  the  Wcibull  Characteristic  life  and  slope 
for  three  different  greases  and  for  windings.  These  tables  list  some  of  the 
significant  test  parameters  such  as  temperature  and  speed.  An  additional 
parameter  that  requires  further  definition  is  indicated  in  the  tables  as 
quality  code.  The  0-Code  refers  to  conrnercial  standards  for  motor  construc- 
tion while  the  1-Code  is  of  military  quality.  Although  considerably  more 
Information  is  available  for  each  data  point,  only  the  more  significant  param- 
eters as  determined  from  analysis  discussed  in  Section  3.2.2  are  tabulated. 


In  order  to  test  the  data  fit  to  a Welbull  distrubitlon,  the  correlation 
coefficient  and  an  F test  was  computed  for  each  test  population  group.  Table 
5 is  a typical  output  from  45  test  ppulations.  Since  a correlation  coefficient 
of  1 is  perfect  correlation  of  data  fit  to  a straight  line,  values  of  0.9  and 
greater  indicate  a high  degree  of  correlation  of  the  data  utilizing  the  Wei- 
bull  Distribution  Function. 

The  correlation  coefficient  is  not  sufficient  for  testing  the  "goodness 
of  fit"  when  the  sample  size  Is  small.  As  the  number  of  test  points  is  re- 
duced, this  measure  of  fit  is  biased  toward  an  unusually  high  value.  In  fact, 
(R)  is  exactly  1 wlien  there  are  two  data  points  being  fit.  An  unbiased  mea> 
sure  of  the  fit  is  provided  by  the  F-test  when  the  sample  size  Is  small.  The 
F -test,  which  accounts  for  the  number  of  data  points  being  fit,  was  used 
in  measuring  the  "goodness  of  fit"  when  the  number  of  test  points  was  stoall. 

By  inspection  the  F-test  (t-test  for  this  single  parameter  fit)  shown  in  Table 
5 reveals  the  Welbull  function  is  a highly  reliable  explainer  of  the  data 
contained  in  forty- five  tests  shown.  Note  that  the  F-value  is  high  (typically 
10  or  more)  for  the  majority  of  tests.  The  significance  level  of  the  Welbull 
fit  exceeds  the  95th  percentile  for  most  of  the  data  displayed. 


The  motor  failure  prediction  models  finally  developed  are  based  on  the 
suanerlzed  data  of  Tables  I through  4.  This  data,  therefore,  was  stored  In 
a digital  computer  and  subjected  to  the  analysis  described  In  Appendix  A and 
sunaarlzed  In  Section  3.2.2  of  this  report. 
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60 

4 

4 

0 

4.62 

6558 

125 

160 

6.35 

60 

3 

3 

0 

10.51 

6044 

125 

1A5 

6.35 

no 

9 

9 

0 

2. 11 

4633 

125 

1A5 

6.  35 

no 

5 

5 

0 

2.96 

4786 

125 

165 

6.35 

55 

10 

9 

1 

5.65 

3687 

72 

112 

6.35 

220 

5 

3 

1 

1.63 

13533 

100 

no 

6.35 

103 

6 

6 

1 

5.33 

27191 

125 

165 

6.35 

55 

6 

6 

1 

2.34 

4577 

50 

65 

6.35 

33 

15 

8 

1 

1.96 

95907 

100 

120 

6.35 

35 

14 

12 

1 

1.45 

40842 

100 

135 

3.18 

32 

30 

20 

1 

3.07 

43975 

100 

135 

3.18 

32 

12 

7 

1 

2.95 

38353 

72 

109 

4.76 

225 

6 

6 

1 

2.94 

42863 

100 

137 

4.76 

225 

8 

7 

1 

4.85 

11646 

125 

162 

4.76 

225 

10 

10 

1 

2.90 

6273 

72 

105 

4.76 

31 

23 

23 

1 

2.05 

24848 
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TEST 

.'O^IP 

BORE 

SPEED 

NUMBER  OF 

GUAL 

WE I BULL 

CHARACTER 

AMB 

BRG 

DIA 

HNDRS 

TEST 

UNITS 

CODE 

SLOPE 

LIFE 

C 

C 

MM 

RPM 

TOTL 

#FAIL 

« 

BETA 

HOURS 

72 

105 

4.76 

31 

12 

11 

1 

2.25 

25136 

72 

105 

4.76 

31 

6 

5 

1 

4.54 

29433 

72 

105 

4.76 

31 

8 

8 

1 

0.80 

32359 

125 

151 

3. 18 

225 

4 

4 

1 

1.33 

4544 

100 

135 

3. 18 

32 

12 

10 

1 

2.18 

41217 

125 

162 

3. 18 

200 

6 

4 

1 

0.87 

7342 

50 

88 

4.76 

200 

24 

4 

1 

0.70 

765796 

72 

117 

3. 18 

205 

4 

4 

1 

8.56 

14307 

100 

135 

3. 18 

32 

12 

8 

1 

1.65 

47431 

100 

^35 

3. 18 

32 

15 

15 

1 

1.94 

29864 

72 

107 

3. 18 

32 

12 

7 

1 

2.08 

39831 

TABL-e:  s 

l__  I !=■£  TErST  Y 

GFtE/svBE-  S 


TEST 

TEMP 

BORE 

SPEED 

NUMBER  OF 

QUAL 

ImEIBULL 

CHARACTER 

AMB 

BRG 

DIA 

HNDRS 

TEST 

UNITS 

CODE 

SLOPE 

LIFE 

C 

C 

MM 

RPM 

TOTL 

#FAIL 

# 

BETA 

HOURS 

125 

157 

3.  IS 

205 

5 

5 

0 

A. 81 

33A1 

125 

167 

3.18 

205 

5 

5 

0 

8.  lA 

A036 

72 

110 

A.  76 

33 

5 

5 

1 

1.72 

A1850 

72 

32 

6.35 

35 

7 

7 

0 

0.98 

17B60 

72 

llA 

6.35 

112 

3 

3 

0 

1.62 

2752 

72 

89 

6.35 

33 

3 

3 

1 

6.00 

26750 

72 

89 

6.35 

33 

3 

3 

1 

3.8A 

8133 

72 

92 

8.10 

33 

5 

5 

1 

5.3A 

1199A 

72 

105 

A.  76 

31 

7 

6 

1 

A.  93 

29966 

72 

105 

A.  76 

31 

7 

7 

1 

A.  58 

36531 

72 

105 

A.  76 

31 

7 

7 

1 

A.  55 

379AA 

72 

100 

3. 18 

105 

A 

A 

1 

1.37 

22081 

72 

110 

A.  76 

33 

6 

3 

0 

1.28 

16938 

72 

llA 

3. 18 

33 

5 

A 

1 

2.18 

173AA 

72 

110 

A.  76 

225 

A 

A 

1 

2.  87 

21322 

72 

110 

A.  76 

33 

6 

3 

0 

1.28 

16938 

72 

92 

8.10 

33 

6 

6 

0 

2.A5 

3936 

11 
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l_  T R’Cf£  "T"EIS3  T 53t-l»»fM/aiFt  V 

GFtEErAoe;  3 

TEST 

TEMP 

B^E 

SPEED 

NUMBER  OF 

QUAL 

WE  I BULL. 

CHARACTER. 

AMB 

BRG 

DIA 

HNDRS 

TEST 

UNITS 

CIX)E 

SLOPE 

LIFE 

C 

C 

MM 

PPM 

TOTL 

#FAIL 

« 

BETA 

HOLIRB 

125 

157 

3.  18 

205 

14 

9 

1 

3.63 

9352 

125 

1E7 

3.18 

205 

32 

30 

0 

2.41 

4054 

125 

177 

3.  18 

205 

37 

28 

0 

3.41 

2854 

125 

16.7 

3.18 

205 

33 

26 

0 

2.61 

3109 

125 

177 

3.  18 

205 

13 

8 

0 

3.91 

2297 

125 

177 

3.  IS 

205 

5 

3 

0 

1.57 

784 

125 

167 

3.18 

205 

17 

8 

0 

2.81 

2799 

125 

167 

3.  18 

205 

10 

4 

0 

3.62 

4338 

125 

177 

3.18 

205 

5 

4 

1 

35.61 

1858 

125 

177 

3.  18 

205 

7 

5 

1 

4.34 

2694 

125 

177 

3.  18 

205 

6 

4 

1 

2.49 

1687 

125 

160 

4.76 

135 

3 

3 

0 

3.85 

1440 

125 

163 

4.76 

225 

5 

5 

0 

1.26 

1588 

125 

163 

4.76 

225 

5 

5 

0 

4.07 

1066 

125 

163 

4.76 

225 

8 

7 

0 

1.24 

537 

125 

163 

4.76 

225 

3 

3 

0 

1.17 

872 

125 

158 

4.76 

110 

6 

6 

0 

1.79 

3740 

125 

163 

4.76 

225 

4 

4 

0 

3.60 

1081 

125 

163 

4.76 

225 

5 

5 

1 

3.51 

2647 

125 

163 

4.76 

225 

4 

4 

0 

0.76 

861 

125 

163 

4.76 

225 

5 

5 

0 

1.26 

1588 

125 

163 

4.76 

225 

5 

5 

0 

4.07 

1066 

125 

163 

4.76 

225 

7 

6 

0 

0.93 

918 

72 

152 

4.76 

165 

4 

4 

1 

0.84 

3237 

125 

163 

4.76 

225 

4 

4 

0 

5.80 

676 

110 

148 

4.76 

33 

3 

3 

0 

9.02 

8669 

125 

150 

6.35 

no 

4 

4 

0 

1.46 

380 

125 

160 

6.35 

60 

4 

4 

0 

1.97 

4345 

125 

214 

6.35 

100 

b 

5 

0 

2.26 

618 

135 

170 

6.35 

55 

5 

5 

0 

3.21 

1130 

135 

170 

6.35 

55 

5 

5 

0 

4.19 

1530 

125 

214 

6.35 

100 

5 

5 

0 

2.90 

591 

135 

201 

6.35 

55 

5 

5 

0 

3.24 

1611 

125 

214 

6.35 

100 

3 

3 

0 

2.37 

1170 

135 

203 

6.35 

60 

5 

5 

0 

3.91 

1868 

125 

155 

6.35 

35 

6 

6 

0 

3.01 

4006 

125 

145 

6.35 

75 

6 

6 

0 

3.32 

2681 

125 

160 

6.35 

60 

7 

7 

0 

4.54 

853 

125 

193 

6.35 

118 

3 

3 

0 

4.48 

186 

125 

193 

6.35 

118 

3 

3 

0 

1.23 

1371 

Ml; 
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VEST 

TEMP 

BmE 

SPEED 

NUMBER  OF 

QUAL 

WEIBU.L 

CHARACTER 

amb 

BHG 

DIA 

HNORS 

TEST 

UNITS 

CODE 

SLOPE 

LIFE 

r 

w- 

C 

MM 

RPM 

TOTL 

#FAIL 

# 

BETA 

HOURS 

72 

140 

S.35 

UX) 

7 

7 

0 

3.01 

301 

I?';* 

ISO 

fc  ".IC 

«-*  • ••.•'A 

50 

S 

s 

0 

l.SS 

14S 

12S 

214 

S.35 

100 

3 

3 

0 

2.73 

584 

las 

214 

S.  35 

100 

5 

5 

0 

l.SS 

39S 

las 

150 

S.  35 

SO 

4 

4 

0 

2.91 

1585 

130 

155 

S.  35 

so 

15 

lA 

0 

4.39 

295 

las 

ISO 

S.35 

SO 

3 

3 

0 

1.14 

6SS3 

las 

ISO 

S.35 

SO 

5 

5 

0 

1.99 

5456 

las 

ISO 

S.35 

SO 

5 

5 

0 

2.04 

4155 

las 

ISO 

S.35 

SO 

15 

15 

1 

2.98 

1417 

las 

17S 

3.  18 

205 

3 

3 

0 

1.08 

1431 

las 

aoo 

3.  18 

210 

3 

3 

0 

1.01 

879 

las 

aoo 

3.  18 

210 

3 

3 

0 

1.57 

831 

las 

175 

3. 18 

205 

5 

5 

0 

3.68 

1746 

las 

132 

3.  18 

205 

3 

3 

0 

0.73 

1444 

72 

105 

4.7S 

31 

15 

14 

1 

3.97 

18574 

72 

105 

4.7S 

31 

15 

14 

1 

2.55 

24438 

72 

105 

4.7S 

110 

5 

5 

1 

1.25 

10604 

lAO 

182 

3.  18 

205 

IS 

4 

1 

1.84 

855 

las 

177 

3.  18 

205 

7 

7 

0 

1.59 

4023 

125 

177 

3.  18 

225 

4 

3 

0 

1.5S 

4925 

las 

1S3 

4.7S 

225 

33 

25 

0 

1.55 

1883 

125 

175 

3.  18 

205 

8 

S 

0 

1.32 

2474 

125 

163 

4.7S 

225 

12 

7 

0 

1.9S 

1038 

125 

150 

4.7S 

225 

7 

7 

0 

1.49 

3286 

125 

1S3 

4.7S 

225 

4 

4 

0 

3.21 

862 

125 

1S3 

4.7S 

225 

6 

5 

1 

1.35 

2127 

125 

175 

3.  18 

200 

5 

3 

0 

1.S5 

3312 

13 
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Ul  ^ IMO  I rMC;  A I L_LJFtiE:0 


TEST 

TEMP 

BORE 

SPEED 

NUMBER  OF 

OUAL 

WEIBULL 

CHARACTER 

AMB 

BRG 

DIA 

HNDRS 

TEST 

UNITS 

CODE 

SLOPE 

LIFE 

C 

t 

MM 

RPM 

TOTL 

#FAIL 

# 

BETA 

HOURS 

las 

177 

3. 18 

205 

4 

4 

1 

1.47 

452 

lAO 

182 

3. 18 

205 

17 

16 

1 

0.98 

1348 

lAO 

1P2 

3.18 

205 

12 

11 

1 

0.74 

726 

125 

177 

3.18 

205 

7 

6 

1 

2.65 

4569 

125 

177 

3. 18 

205 

f. 

5 

1 

1.36 

4359 

lAO 

152 

3.  18 

205 

e 

6 

1 

2.30 

2689 

140 

1B2 

3. 18 

205 

10 

5 

1 

2.41 

2516 

125 

?69 

3. 18 

no 

30 

17 

1 

1.35 

8124 

125 

177 

3. 18 

205 

25 

7 

1 

1.62 

15373 

125 

157 

3. 18 

205 

29 

29 

1 

2.51 

14527 

12  b 

157 

3. 18 

205 

J4 

5 

1 

4.09 

9568 

125 

157 

3. 18 

205 

33 

3 

0 

4.45 

2306 

125 

167 

3. 18 

205 

5 

4 

0 

1.09 

182 

125 

157 

3.18 

205 

13 

3 

0 

0.77 

4702 

125 

157 

3.18 

205 

17 

7 

0 

0.75 

2742 

125 

157 

3. 18 

205 

10 

6 

0 

1.02 

4073 

125 

177 

3. 18 

205 

20 

6 

0 

3.63 

4720 

125 

157 

3.13 

205 

20 

10 

0 

3.77 

4531 

125 

157 

3.18 

205 

5 

3 

0 

9.29 

6545 

125 

157 

3. 18 

205 

6 

1 

1.91 

7429 

140 

1V4 

3.18 

no 

12 

10 

1 

3. 13 

411 

140 

174 

3.18 

no 

14 

14 

1 

3.97 

403 

125 

157 

3.18 

205 

17 

5 

1 

4.  n 

7930 

125 

177 

3.18 

20S 

30 

16 

1 

1.82 

11270 

50 

87 

3. 18 

200 

8 

3 

1 

0.49 

230883 

24 

55 

3.18 

205 

15 

4 

1 

0.51 

716959 

50 

32 

3.18 

205 

IS 

4 

1 

1.04 

60025 

100 

175 

3.18 

32 

30 

10 

1 

1.14 

£0023 

140 

1 ob 

3.18 

115 

5 

5 

1 

18.22 

113 

125 

170 

3.18 

115 

10 

9 

1 

3.22 

1890 

1£5 

170 

3.18 

115 

12 

12 

1 

1.07 

168 

152 

180 

3.  IS 

145 

6 

6 

1 

3.26 

1532 

140 

158 

3.18 

145 

11 

10 

1 

3.70 

8397 

140 

182 

3.18 

205 

16 

9 

1 

1.91 

1476 

125 

153 

4.75 

225 

33 

4 

0 

0.68 

15334 

50 

87 

3.18 

200 

24 

4 

1 

0,50 

999999 

125 

153 

4.75 

225 

5 

5 

0 

4.13 

564 

u 


NUMBER  OF 

WE  I BULL. 

CHAR. 

CORR 

F TEST 

TEST 

UNITS 

CODE 

SLOPE 

LIFE 

COEF 

VALUE 

TOTL. 

tFAIL 

t»1» 

beta 

HOURS 

#« 

14 

9 

1 

3.E3 

93Ba 

.983 

50A.  1 

3S. 

30 

0 

a.  Ai 

AOBA 

.95A 

586.3 

37 

58 

0 

3.  A1 

assA 

.970 

Aa3.0 

33 

5E 

0 

a. El 

3109 

.903 

106.7 

13 

8 

0 

3.91 

5537 

.975 

105.9 

5 

3 

0 

1.B7 

78A 

.95E 

10.7 

17 

8 

0 

a. 81 

5799 

.9A9 

55.5 

10 

A 

0 

3.Ea 

A338 

.989 

91.  E 

5 

A 

1 

35.  El 

1853 

.901 

8.E 

7 

5 

1 

A,3A 

acsA 

.9A9 

57.6 

E 

A 

1 

a.A9 

1E87 

.88A 

7. a 

3 

3 

0 

3.&B 

lAAO 

.908 

A. 7 

S 

S 

0 

i.as 

1588 

.9A0 

53.0 

5 

5 

0 

A.  07 

1066 

.908 

lA.  1 

8 

7 

0 

i.aA 

537 

.976 

103.  5 

3 

3 

0 

1.17 

875 

.995 

118.1 

6 

E 

0 

1.79 

37  AO 

.EE9 

3.5 

A 

A 

0 

3.60 

1081 

.937 

lA.E 

S 

5 

1 

3.51 

aEA7 

.965 

AO.  9 

A 

A 

0 

0.7E 

861 

.8E7 

E.O 

5 

5 

0 

i.ae 

1588 

.9A0 

53.0 

5 

5 

0 

A.  07 

lOEE 

.908 

lA.  1 

7 

E 

0 

0.93 

918 

.813 

7.8 

A 

A 

1 

0.8A 

3337 

.975 

39.5 

A 

A 

0 

5.80 

E7E 

.891 

7.7 

3 

3 

0 

s.oa 

8EE9 

.999 

535.1 

A 

A 

0 

I.AE 

380 

.885 

A. 5 

A 

A 

0 

1.97 

A3A5 

.997 

365.3 

5 

5 

0 

a. as 

E18 

.985 

8A.8 

5 

5 

0 

3.ai 

1130 

.805 

5.  A 

5 

5 

0 

A.  19 

1B30 

.90A 

13.A 

5 

5 

0 

a. 90 

591 

.885 

10.5 

5 

5 

r> 

w 

3.aA 

lEll 

.966 

Aa.5 

3 

3 

0 

a.  37 

1170 

.953 

9.9 

5 

B 

0 

3.91 

18E8 

.968 

A5.6 

E 

E 

0 

3.01 

AOOE 

.957 

A3. 7 

6 

E 

0 

3.3a 

8E81 

.919 

58.0 

7 

7 

0 

A.  BA 

853 

.955 

55.5 

3 

3 

0 

A.  AS 

18E 

.987 

38.  A 

3 

3 

0 

1.53 

1371 

.987 

6.1 

7 

7 

0 

3.01 

301 

.991 

595.3 

S 

B 

0 

1.99 

lAE 

.9A9 

57.6 

3 

3 

0 

a.  73 

58A 

.970 

16.  A 

5 

B 

0 

1.E9 

396 

.9A9 

57.3 

A 

A 

0 

8.91 

ISSS 

.9A5 

17.0 

15 
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3.2  Bearing  Reliability  Model 

As  previously  discussed,  the  reliability  model  to  be  developed  for 
FHP  motors  Is  dominated  by  tvo  failure  modes  bearings  and  windings.  Wel- 
bull  characteristic  life  and  slope  at  different  operating  temperatures  and 
condiltlons  were  calculated  for  each  failure  mode.  Ihe  approach  followed  was 
to  develop  a mathematical  model  of  the  influence  on  life  of  parameters  such 
as  temperature,  speed,  motor  type,  etc.  for  each  failure  mode.  An  expression 
oould  then  be  developed  to  re-comblne  both  failure  modes  Into  a single  motor 
failure  model.  The  following  section  discusses  the  development  of  a failure 
model  for  the  motor  bearing,  commencing  with  a literature  review,  developing 
the  model  parameters  of  concern  and  summarizing  the  bearing  data  analysis 
conducted. 

3.2.1  Literature  Review  of  Bearings 

In  order  to  assist  In  the  data  analysis,  a literature  review  of  grease 
lubricated  bearings  was  conducted.  The  purpose  of  such  a review  was  to 
Insure  that  all  of  the  significant  bearing  parameters  that  might  influence 
bearing  life  were  eiicoinpassed  in  the  bearing  model.  Furthermore,  since  all 
motor  failure  data  waa  run  under  light  bearing  load  conditions,  It  was  de- 
sirable to  obtain  trends  from  literatu:a  of  load  Influences  on  bearing  life. 

Literature  was  also  reviewed  to  obtain  experience  in  the  operation  of  grease 
bearings  at  low  temperature  since  data  bank  data  in  that  range  is  somewhat 
limited. 

( 

3.2. 1.1  Load  Influences  on  Bearing  Life 

It  is  commonly  accepted  that  a lubricant  film  sust  be  developed  between 
rolling  elements  of  a bearing  in  order  to  insure  long-life  operation.  This 
film  is  termed  an  elastohydrodynasiic  lubrication  (EHL)  film  for  rolling  ele- 
ment bearings, 

1 

A teneralised  elaetohydrodynamie  theory  for  film  clearance  in  rolling  i 

contact  bearings  is  attributed  to  Dowson  (22).  He  expresses  film  thickness  H I 
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(5) 

Each  of  the  terma  In  Equation  (5)  are  deacribad  by  additional  aquations 
which  are  not  described  in  this  discussion  since  the  interest  is  in  load  W, 

A number  of  sources  of  the  EHL  film  prediction  analysis  are  available.  This 
expression«  however,  vaa  selected  since  it  clearly  illustrates  the  influence 
of  load  W on  film  thickness.  The. reciprocal  of  load  is  raised  to  the  0.13 
power  compared  to  tho  other  terop  in  the  aquation.  Since  lubricated  related 
surface  distress  is  associated  with  operation  in  the  absence  of  an  EHL  film, 
ball  bearing  life  in  teims  of  surface  distress  is  a fi’netlon  of  film  thickness. 
Loss  of  grease  film  in  a bearing  due  to  load  would  lead  directly  to  reduced 
life.  Equation  (5)  suggests  that  bearing  load  will  not  significantly  in- 
fluence film  thickne&a  until  loads  approaching  the  maximum  ccpability  of  the 
bearing  are  approached. 

Experiments  conducted  by  Foon  (3)  indicate  that  grease  film  thicknesses 
after  run-in  are  generally  less  than  oil  films.  Ihese  experiments  do  verify 
however,  the  existence  of  a grease  film  in  the  contact  region.  Hie  develop- 
ment of  such  a film  is  undoubtedly  a hydrodynamic  type  of  action.  Therefore, 
it  is  anticipated  that  the  influence  of  load  on  film  for  grease  lubrication  is 
similar  to  that  given  in  Equation  (5)  where  speed  contained  in  the  (U)  term 
is  an  important  parameter  since  it  is  raised  to  the  0.7  power.  In  practice, 
Scarlett  (4)  indicates  few  grease  lubricated  bearings  operate  at  loadings 
that  are  critical  and  load  can  generally  be  Ignored.  It  is  stated  however, 
at  higher  speeds  grease  does  not  act  as  a good  coolant  and  load  Influences 
become  important.  Therefore,  heat  generated  at  higher  speeds  Is  important  in 
reducing  lubricant  viscosity  (u^).  New  Departure  (5)  provides  a nomograph 
that  relates  the  life  reduction  due  to  load  on  grease  lubricated  bearings. 

The  nomograph  was  reduced  to  equatlcv-i  form  for  deep  groove  bearings  utilised 
in  fractional  horsepower  motors  where  the  inner  race  rotates,  as  follows: 
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(6) 


Life  multiplier  M « 1 


where;  N 
W 
SP 


RPM 

Applied  load  (lbs) 

Specific  dynamic  capacity  at  33-1/3  RPM  (lbs) 


To  obtain  an  Insight  into  the  influence  of  load  from  this  equation,  consider 
a typical  FHP  motor  bearing,  i.e.:  9 mm  bore,  llte  specific  dynamic  catalog 

capacity  of  this  bearing  is,  SP  ••  AflS  lbs.  Except  for  gear  or  belt  drive 
motor  applications,  the  typical  bearing  load  is  in  the  range  of  2 pounds. 

From  Equation  (6): 

M 1 IJL\  2 
” " ^ ’ [933  { 485 

The  normal  speed  range  for  a FHP  motor  is  900  RPM  to  20,000  RPM  maximum. 

The  life  multiplier  ranges  from  M • .996  to  .966. 

The  correction  to  life  over  the  operating  range  is  less  than  four  percent 
and  could  be  ignored.  If  however  the  bearing  is  loaded  to  50  percent  of 
specific  load,  life  is  reduced  51  percent  at  900  rpm  and  approaches  zero  at 
2500  rpm.  Therefore,  high  loads  introduce  significant  reductions  in  bearing 
life. 


In  order  to  check  the  accuracy  of  this  load  factor,  a comparison  was 
made  with  the  following  data  taken  from  Reference  6: 


Bearing  204  (20  mm  bore) 

Speed  10,000  RPM 

Amb  232°C 

Axial  Load 

Life 

Predicted  Life 

40  lbs. 

1000  hrs. 

1000  hrs. 

230  lbs. 

100  hrs. 

293  hrs. 

330  lbs. 

10  hrs. 

10  hrs. 
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The  prediction  of  life  assumed  1000  hours  life  at  40  lbs.  load  and  predicted  the 
reduction  in  life  at  230  and  330  lbs.  The  correlation  of  life  reduction  with 
load  checks  quite  well  with  Reference  (6).  Other  direct  correlations  of  load 
Influences  on  life  at  constant  temperature  and  speed  were  not  found  and 
therefore  Equation  (6)  was  selected  as  a suitable  expression  to  relate  load 
Influences  on  life. 


Reference  (1)  developed  a relationship  of  the  influence  of  temperature 
of  the  bearing  on  life  In  the  form: 

Log  life  “ ^ + B (7) 

where:  A and  B are  constants 

r ■ Bearing  absolute  temperature 


It  Is  desirable  to  describe  the  influence  of  load  on  life  In  logarithmic  form 
to  provide  compatibility  with  Equation  (7),  A close  approximation  of  Equation 
(6)  in  log  form  is: 


.OOIN 


r 

w_ 

SP 


1.5 


= Log  M 


(8) 


The  generalized  life  formula  of  Equation  (7)  would  thus  be  modified  by 
Equation  (8)  i.e,: 

1.5  (9) 

Log  life  - £ -r  B - .OOIN  I — 


3.2. 1.2  Low  Tempecatuie  Influences  on  Bearing  Life 

As  the  operating  aidbi‘?int  temperature  of  a bearing  is  decreased,  eventually 
a temperature  must  be  attained  where  the  lubricant  can  no  longer  perform  its 
function  because  it  approaches  a solid  state.  It  would  be  anticipated  that  the 
lubrication  mechanism  does  not  cease  ct  ai.  exact  temperature  but  is  subject  to 
a transition  range  where  lubrication  becomes  marginal.  Following  this  assump- 
tion, a temperature  must  exist  where  the  life  of  the  bearing  no  longer  in- 
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creases  but  begins  to  decrease  with  decreasing  temperature.  The  literature 
/review  provided  no  available  data  of  grease  lubricated  bearing  life  at  low 
temperature.  Low  temperature  testing  by  Crisp  and  Wells  (7)  resulted  in 
occasional  damage  in  the  bearing  after  short  periods  of  operation  at  ambients 
in  the  range  of  -38 °C,  They  concluded  that  damage  noted  in  the  sliding  regions 
of  the  bearing  was  indicative  of  marginal  lubrication;  however,  there  was  no 
indication  of  the  anticipated  life.  Other  investigators  (8,9)  have  been  con- 
cerned with  bearing  torque  at  low  temperatures.  In  all  torque  measurement 
tests  of  greases  at  low  temperature,  a temperature  is  reached  where  torque 
increas<»s  rapidly  with  decreasing  temperature.  If  Crisp  (7)  is  correct  that 
the  increasing  torque  level  may  be  associated  in  part  with  marginal  lubrication, 
then  the  transition  temperature  to  high  torque  operation  may  be  indicative  of 
reduced  lubrication  and  bearing  life.  At  very  low  bearing  speeds  Increasing 
torque  appears  to  be  a function  of  lubricant  properties  (9)  and  lubricant 
shear  resistance.  As  speeds  are  increased,  the  onset  temperature  of  increas- 
ing torque  increases  (7)  and  at  DN  values  of  7000  a transition  temperature  of 
16°C  outer  race  temperature  was  measured.  If  this  is  due  in  part  to  marginal 
lubrication,  then  one  might  expect  bearing  life  of  FHP  motors  to  begin  to  fall 
off  at  bearing  temperatures  in  the  range  of  16^C.  Since  literature  conceimed 
with  life  at  low  temperature  of  grease  packed  bearings  is  so  limited,  consider- 
able conjecture  is  associated  with  estimating  the  lower  temperature  bound  of 
Equation  (7).  Although  actual  life  test  data  from  the  data  bank  is  used  to  es- 
tablish the  useful  range  of  Equation  (7),  it  is  desirable  to  verify  these 
trends  from  available  literature.  The  trends  worthy  of  note  from  this  review 
are  that  reduction  in  bearing  life  with  temperature  may  occur  at  rather  high 
ambient  temperatures,  and  the  transition  point  is  a function  of  beariqg  /^peed 
and  lubricant  characteristics.  The  reader  is  referred  to  Section  3. 2. 2. 4 for 
a discussion  of  experimental  results  regarding  low  temperature  life  of  bear- 
ings. 


3.2.2  Development  of  Bearing  Life  Prediction  Model 

The  bearing  life  prediction  model  was  developed  from  the  data  tabulated 
in  Tables  1 through  3 utilizing  a multivariable  regression  analysis  computer 
program.  This  program  and  its  capabilities  is  described  in  Appendix  A. 
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Reference  I established  bearing  grease  life  as  strongly  Influenced  by 
temperature  In  the  form  of  Equation  (7) : 


Log  life  - f - B 


(7) 


This  equation  linearizes  the  relationship  between  life  and  temperature  per- 
mitting the  use  of  linear  regression  techniques  to  establish  the  coefficients 
A and  B.  To  Illustrate  the  relationship  of  Equation  (7) » Figure  1 Is  a plot 
of  the  log  of  characteristic  life  of  all  data  from  the  data  bank  Including 
the  data  of  Tables  1 through  3 as  a function  of  the  reciprocal  of  temperature. 
Without  Introducing  any  additional  variables,  the  strong  dependency  of  life 
on  temperature  Is  apparent.  A linear  regression  of  this  data  results  In  a 
solution  of  Equation  (7)  as  : 

Log  lifeg  - - 1.908  (10) 

where:  Life  ■ Weibull  characteristic  life  (a_)  in  hours 

^ o 

T • Absolute  temperature  K 


This  equation  explains  54  percent  of  the  data  scatter  of  Figure  1.  A further 
reduction  In  scatter  can  be  effected  by  Including  additional  variables  that 
Influence  life.  This  Implies  that  the  coefficient  B of  Equation  (7)  is  com- 
posed of  a number  of  additional  variables. 

The  variables  entered  Into  regression  included; 


Bearing  operating  temperature 

Bearing  speed 

Bearing  site 

Number  of  motor  poles 


Quality  code  and  grease  type 
Winding  temperature  rise 
Motor  voltage 
Motor  winding  site 


The  standard  F-test  was  used  to  determine  the  significance  of  correlation 
between  the  observed  teat  failure  livea  and  the  various  operating  parameters 
listed.  The  dependent  variafbles  on  bearing  life  were; 
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Bearing  temperature 

Bearing  size  (bore  diameter  In  mm)  times  speed 
Quality  code 
Grease  type 


In  separating  grease  type  In  the  regression  analysis,  it  was  c erved 
that  the  slope  coefficient  A of  Equation  (7)  varied  for  each  grease  ranging 
from  1730  to  2441,  Although  one  might  anticipate  that  different  greases  might 
have  a different  slope  of  life  versus  temperature,  the  distribution  of  fail- 
ures over  the  temperature  range  tor  two  greases  was  skewed.  Tills  tended  to 
reduce  the  reliability  of  the  slope  prediction  for  Individual  grease.  As  a 
rest' t,  the  slope  for  all  greases  was  maintained  at  2342  from  Equation  (10) 
and  the  Influence  of  quality  code,  speed,  and  grease  type  was  evaluated. 


Using  Equation  (4)  for  Illustrative  purposes,  the  equation  Is  rewritten 
in  the  following  form  to  obtain  the  Influence  of  quality  code; 


Log  lifSjj 


qual  code:  0\ 
qual  code:  1 


<11) 


The  mean  and  standard  deviation  of  x.  and  x values  are  calculated.  These 

I 2 

corrections  are  then  transposed  to  the  left-hand  side  and  the  Influence  of 
speed  determined  In  the  same  manner;  l.e. 


Log  lifSg  - 


m 4- 


(12) 


A regression  Is  run  on  this  equation  to  determine  the  coefficients  and 
x^.  These  coefficients  are  then  transposed  to  the  left  side  of  the  equation 
and  the  Influence  of  each  grease  obtained: 


As  the  coefficients  are  determined,  the  predictive  equations  can  be  rewritten: 
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(10) 


Log  llfCg 
Log  Ilf eg 
Log  lifeg 
Log  llfeg 


2342 

"b 


2342 

"b 


2342 

^B 


2342 

T 

B 


1.91 

1.91  + q 

1.91  + q - eDN 

eDN  + q - k 

8 


q ■ quality  code 


e " constant: 


K • grease  constant 
& 


(14) 

(15) 


(16) 


The  data  «as  tested  against  each  of  these  equations  and  percentage  of  actual 
lives  falling  withlr  one  standard  deviation  were  calculated.  A standard 
deviation  in  log  form  of  .35  for  the  intercept  coefficient  B of  Equation  (7) 
was  used.  Table  6 summarizes  the  equations  and  data  fit. 


The  general  equation  for  bearing  characteristic  life  (hrs)  may  be  written 
In  the  form: 


Log  11  fag  ■ 


+ q - 4.32  DN  X 10“®  - K 


g 


where:  q 


DN 

K 

8 

c 


quality  code 

comaercial  quality  q • -.27 
Military  spec,  quality  q ■ +.12 
Bearing  bore  (mm)  x speed  (rpm) 
Greaae  constant  per  Table  7 
4.32  X 10  ^ a constant 

TABLE  7 

Greaae  Constant  K 

8 


Grease 

1 

2 

3 


K 

8 


(17) 
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1.35 

1.55 

1.74 


TABLE  6 


Test  of  Bearing  Failure  Data 
Versus 

Predictive  Equations  for  One  Standard  Deviation** 


Percent  of  Data  Within 
Equation  One  Standard  Deviation 


A. 

2342 

^B 

1.91 

69.85 

B. 

2342  I 

[2.03  LQG*) 

70.59 

"b  ■' 

[l.64  UQG  1 

C. 

2342 

[I'.M  uq^)- 

74.26 

C 

Ito 

hIw 

00 

Ifo 

1 

[ i as  wj-  ‘-32  « w'*™ 

'l.50 

-'1.70 

\1.89 

GRl' 

GR2 

GR3/ 

75.00 

LQG  • low  quality  grade 
UQG  ■ upper  quality  grade 

One  standard  deviation  was  found  to  be 
.35  on  the  log  life  (hrs)  scale 


The  data  bank  did  not  contain  sufficient  bearing  load  variation  to  in- 
clude load  influences  on  life  predictions.  It  was  necessary  to  develop  a 
load-life  modifier  from  available  Information  found  in  literature.  As  dis- 
cussed in  Section  T.2.1.1,  Equation  (4)  and  (5),  Equation  (17)  may  be  rewritten 
to  include  load  effects. 


Log  Ufe^- 


+ q - 4.23DN  x lO"®  - Kg  - .OOIN 


W 

SP 


ll.5 


(18) 


The  useful  temperature  range  of  this  equation  varies  from  the  maximum 
recommended  temperature  by  the  grease  manufacturer  to  a temperature  in  the 
range  of  IS^C.  Section  3. 2.2.4  should  be  referred  to  for  development  of  the 
low  temperature  life  characteristics. 

3. 2. 2.1  Bearing  Temperature  Rise 

Equation  (18)  was  developed  for  predicting  life  requiring  considerable 
knowledge  regarding  the  bearing;  l.e.  the  type  of  grease,  the  bearing  temper- 
ature, bearing  size  and  speed.  Since  bearing  temperature  Is  the  dominant  life 
Influence,  accurate  determination  of  this  temperature  is  quite  Important. 

Motor  bearing  temperature  is  a difficult  parameter  to  ascertain.  Wind- 
ing temperature  rise  above  and>ient  temperature  is  often  a good  indicator  of 
bearing  temperature.  Generally,  bearing  temperatures  (outer  race)  run  lO^C 
cooler  than  windings.  In  some  type  of  laotors  however,  the  bearings  may  run 
hotter  than  the  winding  depending  on  their  cooling  and  location  in  the  motor. 
Table  8 Illustrates  soma  typical  maasurements  oiade  on  a number  of  different 
type  motors.  As  noted,  bearing  temperature  rise  above  ambient  ranges  from  15*^ 
to  ISO^C.  The  distribution  of  outor  winding  temperature  rise  taken  by  resist- 
ance swasurements  for  a hundred  and  seventy-five  different  fan-cooled  FHP 
motors,  ranging  from  polyphase  to  capacitor-run  and  shaded  poles  is  pre- 
aentad  In  Figure  2.  The  mean  value  for  winding  rise  la  36.4°C.  The 
mean  value  for  bearing  temperature  rise  would  be  26.4°C.  The  standard  devia- 
tion is  approximately  16*’c.  Therefore,  one  can  anticipate  that  a 68  per- 
cent probability  exists  that  the  bearing  temperature  will  fall  in  the  range 


STEP  f'lORWALIZED  DATA  HISTt^RAM 

A. 000  X 

8.000  X 

12.000  X 

16.000  XXXXXXXXXXXX 

20.000  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

2A.000  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

28.000  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

32. 000  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

36.000  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

40.000  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

44.000  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

48.000  xxxxxxxxxxxxxxxxxxxxxxxxx 

B2. 000  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

56.000  XXXXXXXXXXXXXXX 

60.000  xxxxxxx 

64.000  xxxxxxx 

68.000  XXXXX 

72.000  XXXXXXXXXX 

76.000  xxxxxxx 

80.000  XX 

84.000  XX 

88.000  XX 

92.000  X 

96.000  XX 

00.000  AREA  = 175.0  STD  DEV  » 015.8507  PEAK  » 020. 


Figure  2 

Motor  Ttaiperature  Rise  Dittrlbutlon 
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MP 


of  10°  to  43°C.  Recognising  that  there  is  approximately  a 50  percent  reduc" 
tlon  In  life  for  every  10°C  rise  in  temperature,  this  range  la  quite  broad. 

The  conaervatlve  approach  when  the  temperature  is  not  known  is  to  use  a bearing 
temperature  rise  of  40°C.  This  temperature  agrees  with  the  data  failure  bank 
analysis  which  Indicated  a mean  bearing  temperature  rise  of  37.9°C.  It  must 
be  recognized  for  high  slip  notors  and  special  application  motors  this  temper- 
ature stay  not  be  conservative. 


TABLE  8 

Typical  Motor  Temperature  Rise  Date 


Winding 

Bearing 

T 

T 

Shaded  Pole  Motor 

57®C 

42°C  (Pull  Load) 

Capacitor  Kun  Motor 

160°C 

150°C  (Full  Load) 

(Totally  Enclosed 

No  Cooling  Fan) 

77°C 

66°C  (I.ight  Load) 

Capacitor  Run  Motor 

95°C 

82 °C  (Pull  Load 

(Vented  — No  Cool- 

ing Fan) 

Hysteresis  Motor 

70°C 

46°C  (Pull  Load) 

Polyphase  (fan  cooled) 

46°C 

23°C  (Light  Load) 

Polyphase  Motor 

74°C 

65°C  (No  Load) 

(Totally  Enclosed  -- 

No  Cooling) 

Inside-Out  Polyphase 

32°C 

50°C 

(Fan  Cooled) 

High  Slip  Motors 

30°C 

15°C  (Low  Slip  (5 

(Fan  Cooled) 

50°C 

90°C  (80%  Slip) 

Polyphase  Motor 

(Fan  Cooled) 

20°C 

15°C 

3. 2. 2. 2 Ambient  Temperature  Influences  on  Bearing  Life 

It  is  often  difficult  for  the  user  of  FHP  notors  to  obtain  informtion 
regarding  the  bearing  temperature  rise  over  at^ient  temperature.  It  would  be 
desirable  therefore,  to  simplify  Equation  (18)  and  relate  bearing  life  in 
terms  of  asblent  trapOvature.  In  addition,  the  improveiiMKit  in  fit  of  data 
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within  a standard  deviation  by  correcting  for  speed,  quality  and  grease  from 
Table  6 eas  only  five  percent,  therefore,  the  simple  model  of  Equation  (10) 
is  recommended  vhen  detailed  knowledge  of  the  motor  Is  not  available. 

As  a further  simplification,  in  lieu  of  entering  bearing  temperature  Into 
regression,  motor  ambient  temperature  was  utilised  since  this  parameter  was 
statistically  significant.  Equation  (7)  could  be  rewritten  using  ambient 
temperature  as : 

Log  life_  - . 2.534  (19) 

When  the  data  in  the  data  bank  was  tested  against  this  model,  63.97  of 
the  data  fell  %fithln  one  standard  deviation.  This  compares  with  69.85  percent 
using  bearing  teuperaturo.  Due  to  the  possible  error  introduced  in  estimating 
bearing  temperature.  Equation  (19)  was  selected  as  the  -general  simplified  siodel 
for  estlmeting  bearing  life.  Under  normal  load-speed  conditions  found  in  FHP 
motors,  this  equation  will  yield  reasonable  estimates  oi  life.  For  unusual 
load  -speed  on  motor  designs,  the  more  comprehensive  Equation  (18)  should  be 
used. 

3.2. 2. 3 Estimating  Bearing  Sire  and  Speed 

Vlhen  utilising  Equation  (18)  for  predicting  bearing  life.  It  Is  often 
possible  to  estimate  bearing  sise  from  the  owtor  Input  power  level.  This  in- 
fotmacion  can  be  utilised  to  calculate  the  DN  value  of  Equation  (18).  Table 
9 Is  a sumsuiry  of  average  bearing  else  after  reviewing  over  two  hundred  dif- 
ferent motors  of  various  horse-power  level.  Bearing  bore  dlasMter  Is  tabulated 
as  a function  of  motor  Input  power. 


TABLE  9 


Baarlna  Sl««  Vs.  Motor  Input  Poocr 


Bearing  Bora  (D) 
(nn) 


Power  Input 
(watts) 

f » 60  Hs ff  ■ 400  tte 


20 

250-2500 

13 

100-  500 

200-1000 

8 

20-  200 

150-  400 

6.35 

10-  50 

80-  130 

4.76 

7-  12 

60-  90 

3.18 

5-  8 

20-  45 

Motor  speed  can  be  estimated  from 


N(RPM)  - 


(20) 


where:  f 

P 


line  frequency,  Hz 
number  of  poles 


3. 2. 2. 4 Bearing  Failure  Data  Analysis  for  Low  Temperaturo 


In  order  to  establish  the  low  temperature  characteristics  of  the  bear- 
ings, the  cyclic  tests  of  Table  10  were  utilized  to  estimate  life.  Thermal 
cycles  for  all  tests  were  not  Identical  and  the  duty  cycles  are  sumarlzed 
In  Table  11. 
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TEST  TEMP 

BORE 

SPEED 

TYPE 

DIA 

HNDRS 

**♦ 

MM 

RPM 

AT 

3.  IB 

205 

TT 

A.7e 

110 

TT 

4.76 

110 

AT 

4.76 

205 

TC 

4.76 

225 

AT 

3. 18 

200 

AT 

3. 18 

200 

TC 

6.35 

90 

TC 

6.  35 

90 

TC 

6.35 

90 

TC 

6.35 

90 

AT 

6.35 

105 

AT 

6.35 

105 

AT 

6. 35 

110 

AT 

6.35 

60 

TC 

8.10 

106 

AT 

3.18 

205 

AT 

3.18 

205 

NUMBER  OF 
TEST  UNITS 

QUAL 

CODE 

NEIBULL 

SLOPE 

CHARACTER 

LIFE 

TOTL  RFAIL 

« 

BETA 

HOURS 

40 

29 

1 

2.66 

5533 

, 5 

5 

1 

0.71 

9328 

5 

5 

1 

1.25 

10604 

5 

3 

1 

0.94 

4822 

5 

4 

1 

2.36 

6653 

30 

3 

1 

1.37 

17094 

10 

8 

1 

2.79 

13931 

6 

6 

1 

:..80 

1758 

4 

4 

1 

2.09 

771 

3 

3 

1 

6.03 

1879 

3 

3 

1 

2.80 

1516 

5 

5 

1 

1.10 

1167 

5 

5 

1 

1.16 

2854 

10 

10 

1 

2.89 

2738 

8 

8 

1 

0.73 

1902 

3 

3 

1 

2.32 

13965 

13 

4 

1 

3.28 

11966 

11 

9 

1 

2.73 

19518 

X^BLE  11 

Thermal  Cycling  Testg 
Duty  Cvclee 


Teet  Numher 

Thermal  Cycle 

1106 

8 hrs.  off 

4 hre.  at  -42*^C  to  +55°  at  3°  per  min. 

656,  622 

39  hrs.  at  -54°C 

241  hrs.  at  +93°C 
remaining  at  +22°C 

1384,  1324,  1456 

1377,  1374 

2 hrs.  off 

2 hrs.  at  -55°C  to  +100°C 

1492 

150  hrs.  at  +25°C 

50  hrs.  at  “55°C 

100  h!M.  at  +125°C 
200  hrs.  at  +75°C 

700  713,  728,  313 

100  hrs.  at  +25°C 

50  hrs.  at  -5j°C 

250  hrs.  at  +55°C 

100  hrs.  at  +85°C 

1533,  1544 

2 hrs.  off 

1 hr.  at  -65°C 

1.5  hrs.  at  -65°C  to 
1.5  hrs  at  71°C 

71°C 

1556,  1517 

1.5  hrs.  off 

2.5  hrs.  at  -55°C 

1.5  hrs.  at  -55°C  to 
1.0  hr.  at  +72°C 

+72°C  ' 

765 

15  hrs.  at  +72°C 

1 hr.  at  +72®C  to  23°C 
7 hre.  «t  23®C 


32 


From  the  characteristic  li£e  o£  Table  10  and  the  duty  cycles  of  Table  11, 


the  nunber  of  thermal  cycles  prior  to  failure  could  be  determined.  The  life 
at  bearing  temperatures  of  50^C  or  greater  could  be  computed  using  Equation 
(18)  previously  developed.  The  characteristic  life  at  lov  temperature  could 
then  be  solved  from  Equation  (21) . 


^hrs.  at  T, 


hrs.  at  T„ 


No 


. of  cycles  (N)  .t  T^  Life  at  T2 


(21) 


TABLE  12 


Lov  Temperature  Bearing  Life 

Teat  T (°C)  Characteristic 

NSber ^ Life  (Hrs) 


1106 

+10 

304 

656 

-21 

46 

622 

-21 

50 

1384 

+28.5 

2041 

1496 

-17 

1112 

1324 

+22.5 

5432 

1456 

+22.5 

4203 

700 

-35 

382 

713 

-35 

80 

728 

-35 

406 

913 

-35 

323 

1533 

-45 

220 

1544 

-45 

635 

1556 

+10 

833 

1517 

+10 

574 

765 

+20 

4628 

1377 

+28.5 

4225 

1374 

+28.5 

9315 

The  resulting  low  teapereture  characteristic  life  determined  from  this 
equation  Is  tabulated  In  Table  12.  The  points  are  also  plotted  on  log- 
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reciprocal  temperature  paper  in  Figure  3.  The  solid  line  is  a fit  of  data 
between  lO^C  and  30°C.  The  equation  of  the  line  is: 


Log  lifftg 


19.7 


(22) 


Considerable  scatter  in  the  data  exists  below  10*^0  and  it  was  decided 
to  bound  the  minimum  life  at  300  hours.  Life  at  low  temperatures  therefore, 
could  be  expressed  in  the  form 


Characteristic  Hfeg  "10  + 


300  hrs. 


(22) 


This  equation  is  shown  as  the  dotted  line  in  Figure  3. 


From  previous  discussions  (Section  3.2. 1.2)  regarding  low  temperature 
life,  it  is  anticipated  that  the  type  grease  and  speed  will  have  considerable 
Influence  on  the  prediction  model  at  low  temperature.  It  is  suspected  that  the 
major  influence  is  to  shift  the  solid  line  of  Figure  3 to  the  left  or  right  and 
the  slope  of  the  line  will  not  be  significantly  changed. 


Since  an  expression  was  derived  to  predict  life  of  the  bearing  as  a 
function  of  ambient  temperature  in  the  higher  temperature  region,  it  is  de- 
sirable to  obtain  a similar  expression  for  low  temperature  life  as  a function 
of  ambient  temperature.  Utilizing  ambient  temperature,  the  equation  derived 
for  life  was: 


log  llf«3" 


-4500 

‘^amb 


+ 20 


(23) 


Setting  the  lower  life  bound  of  300  hours  Equation  (23)  was  modified  to: 


T. 


+ 20! 


tlfij  - lO' 


* + 300  hours 


(24) 
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Combining  the  low  temperature  and  high  temperature  models,  Figure  3 can 
be  modified  to  Include  the  bearing  life  model  over  the  entire  temperature  range. 
The  combined  life  is  illustrated  in  Figure  4. 

It  is  not  expected  that  Che  transition  of  Ufa  from  an  increasing  to  a 
decreasing  function  is  as  abrupt  as  thcj  solid  lines  in  Figure  4.  Therefore,  the 
equations  for  high  temperature  life  and  low  temperature  life  were  combined  in  the 
fonn  of  Equation  (25)  and  shown  as  the  solid  and  dotted  lines  of  Figure  4. 

1 . . A-.  . 

Life  Life  (High  Temp)  Life  (Low  Temp)  (25'' 

Two  models  have  been  developed  for  predicting  bearing  life.  The  more 
comprehensive.  Equation  (18'> , requires  considerable  knowledge  regarding  the 
motor.  Ihe  second.  Equation  (19),  has  a prediction  estimate  accuracy  11  per- 
cent lower  than  Equation  (18);  its  use  requires  knowledge  of  the  ambient 
tenperature  only.  These  two  equations  can  be  modified  to  include  low  temper- 
ature influences  using  Equation  (25)  with  the  follotTlng  overall  bearing  model: 


Life 

- 2.537 

10 

^ amb 

I ^amb 

300  + 10 

L « 

(27) 


The  resulting  life  is  the  Weibull  characteristic  Ufa  or  63,2  percentile 
point  of  the  distribution.  Itm  probability  is  that  63.2  percent  of  the  popu- 
lation will  have  failed  at  this  life. 
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3. 2. 2. 5 W^.J.clull  Slope  for  Bearing  Life 

Equations  (26)  and  (27)  predict  tjpe  bearing  Welbull  characteristic  life 
(ogX  In  order  to  determine  the  life  at  any  other  percentile,  It  Is  neces- 
sary to  establish  the  Welbull  Slope  (Pg)*  In  reliability  prediction  for 
electronic  components,  the  mean  life  Is  of  Interest  In  order  to  determine 
system  failure  rates  for  constant  failure  rate  components.  Although  bearings 
are  an  Increasing  failure  rate  component,  a reasonable  eatlmate  of  the  mean 
life  Is  the  50  percentile  point  of  the  Welbull  distribution.  The  average 
failure  rate  at  this  percentile  Is  suggested  as  a means  of  comparing  bearing 
failure  rate  with  constant  failure  rate  components. 


In  addition  to  establishing  bearing  failure  rates.  It  Is  also  desirable 
on  occasion  to  determine  the  10  percentile  point  often  referred  to  as  the 
life  In  bearing  analysis.  In  order  to  determine  life  at  any  percentile. 
Equation  (18)  or  (19)  can  be  modified  with  knowledge  of  the  Welbull  slope 
(Pj^)  as  follows: 


Log  life  + 


B 


Log  In 


l-F(t) 


I + B 

T 


(28) 


where  E(t)  • failure  percentile 


The  distribution  and  mean  of  the  Welbull  slope  from  Tables  1,  2,  and  3 
were  used  to  establish  the  average  slope  for  use  In  Eqisation  (28).  Figure  5 
is  a histogram  of  the  distribution  of  the  Welbull  slope  of  bearing  failures. 
The  mean  of  the  dlstriburlon  is  2.8789  with  a standard  deviation  of  1.99.  The 
naan  value  of  2.878  was  selected  as  Welbull  slope  P for  bearing  failures. 


3. 2. 2. 6 Bearing  Life  — Grease  Constants 

In  addition  to  the  three  basic  greases,  analysis  was  conducted  on  21 
additional  greases  to  establish  the  grease  constant  K^.  The  values  of  grease 
conetant  ranged  from  1.37  to  2.35  with  a naan  value  of  1.825.  A nanber  of  these 
constants  ware  developed  from  a small  test  population  and  therefore  a question 
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•rises  «s  to  the  rellablltcy  of  the  estimstes.  The  population  o£  test  units 
for  three  greases  was  sufficiently  large  to  Include  reporting  the  values  of 
the  constants. 


The  grease  constants  vere  developed  from  Equation  (18);  l.e. 
+ 27878  Umri  ■ ^ “•  * 


- ■'g  - •ooi"  [fe] 


1.5 


The  resulting  K 


g 


values  for  all  greases  can  be  sunanarlsad  as  follows: 


Grease 

Oil 

Ihlckener 

Mil  Snec. 

Max.  Temo. 

K. 

K . 

1 

dies ter 

Sodium  and  solid 
Lubricant 

MrL-G-3278A 

170 

1.35 

2 

dies ter 

Lithium 

MIL-G-3278A 

120 

1.55 

- 

Silicone 

Lithium 

150 

1.74 

4 

Mineral 

Sodium 

MIL-G-18709A 

150 

1.41 

5 

Silicone 

Lithium 

MIL-L-15719A 

177 

1.81 

6 

Syn. 

H^rocarbon 

Non- soap 

MIL-G-81322 

170 

1.74 

3.*’  i-:.. 

■ (.or  I's  ng  Reliability  Mbdel 

> 

Winding  failures  are  the  second  dominant  failure  mode  of  FXP  motors. 
Table  4 •ummsr''«e  winding  failures  as  a function  of  teoparatura.  By  con- 
ducting multi  . .able  regression  analysis  using  variables  associated  with 
windings,  a reliability  model  was  developed  for  prediction  of  winding  life. 
Tfia  following  sections  discuss  tba  literature  survey  conducted  on  windings  as 
wall  as  davelo^Mnt  of  the  predictive  model. 
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3.3.1  Literature  Review  of  Electrical  Failures 


There  la  general  agreement  between  investlgatora  of  electrical  Insulation 
life  that  chemical  changes  of  the  Insulation  are  dominant  factors  In  Insulation 
deterioration.  It  Is  well  known  that  chemical  raactlon  rates  vary  with  the 
temperature  according  to  the  Arrhenius  Equation  (29) : 

In.u>..tlon  Ilf.  - - A (29) 


or  Log  life 


log  A + I 


This  Is  the  general  form  that  Is  used  to  express  insulation  life.  In  review- 
ing literature  data,  life  was  transposed  into  the  form  of  Equation  (29). 

Most  testing  done  on  Insulation  systems  are  not  conducted  directly  In  the  motor 
but  are  run  In  test  units  called  oiotorettes.  To  accelerate  Insulation  degrada- 
tion, tests  usually  Include  cycles  of  100  percent  humidity  and  vibration.  Data 
presented  In  Reference  10  Indicate  that  the  humidity  cycling  Introduces  between 
4 and  5 to  1 reduction  In  life.  Beebe  (11)  suggests  that  actual  winding  life 
In  motors  exceetbmotorette  life  by  4 to  1.  This  Implies  that  humidity  Is  the 
doodnant  cause  of  life-reduction  In  motorette  tests. 

Table  13  suimiarlses  the  test  data  and  results  reported  in  the  literature: 

moi 


Utarature  Review  of  Temperature  Influence  on  Life 


Halbull 

Slooe 

Insulation 

Clase 

Tyne  Unit 

Husddlty 

Included 

Reference 

^ - 10.7 

T 

m 

B 

Hotoratte 

No 

12 

- lO.M 

- 

R 

Motoretta 

Ho 

12 

^ - 4.615 
T 

«» 

A 

Motor 

No 

ly 

(continued) 
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T«bl«  13  (continued  from 

Wolbull 
Ava.  Life Slop# 


mi 

T 

- 6.95 

- 

mi 

T 

• 8.62 

- 

7474 

T 

- 11.97 

m 

7583 

T 

- 12.2 

- 

im 

T 

- 10.0 

- 

mjL 

T 

- 11.92 

- 

mi 

T 

- 7.02 

•ii 

mi 

T 

- 5.03 

3.0 

im 

T 

- 10.0 

1.6 

im 

T 

- 5.24 

- 

T 

- 8.9 

as 

t 

T 

- 7.02 

- 

T 

- 6.06 

- 

mu 

T 

- 4.96 

so 

T 

- 6.17 

4.5 

T 

- 7.398 

- 

prevloui  page) 

Insulation 

Class  Type  Unit 
A Motorette 

A Motorette 

H Motorette 

H Motorette 

H Motor 

A Motorette 

B Motorette 

A Motor 

B Motor 

A Motorette 

A Motorette 

< A Motor 

A Motor 

A Motor 

A Motorette 

A Motorette 

(continued) 


Humidity 

Included  Reference 


No 

13 

Yes 

13 

Yes 

13 

Yes 

13 

Yes 

14 

Yes 

14 

Yes 

14 

Yes 

15 

Yes 

15 

Yes 

10 

No 

16 

Yes 

16 

Yes 

17 

Yes 

17 

Yes 

18 

Yea 

18 
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Table  13  (continued  from  previous  pag«) 


Hm 

ill 

Welbull 

Slope 

Insulation 

Class 

ipumii 

Humidity 

Included 

Reference 

4406 

T 

- 7.415 

6.8 

A 

Motorette 

Yes 

18 

mi 

T 

- 8.996 

- 

A 

Motorette 

Yes 

11 

4551 

T 

- 7.21 

- 

B 

Motorette 

Yes 

11 

6954 

T 

- 12.17 

- 

F 

Motorette 

Yes 

11 

4994 

T 

- 8.71 

. \ 

A 

Motorette 

Yes 

19 

4785 

T 

- 8.575 

- 

A 

Motorette 

Yes 

19 

1123 

- 8.79 

A 

Motor 

Sc 

20 

T 


A considerable  spread  exists  In  literature  regarding  the  experimental  life 
of  windings.  Ihe  I'eported  reasons  for  the  variation  Is  due  to  difficulty  in 
controlling  test  cycles  such  as  huiiildlty«  variations  in  the  Ittsulatlng  materials 
tested  In  each  winding  class,  different  voltages  applied  to  windings  and  varla> 
tlons  In  interpretation  of  failure  when  a failure  occurs. 

One  further  concern  noted  In  the  available  life  data  Is  the  temperature 
tenge  tued  to  establish  life  characteristics.  The  NEMA  specified  winding 
classes  A,  B,  F and  H correspond  In  general  to  IBBB  classes  lOS^’c,  130**C, 

ISO^C  and  180°C.  These  tas^Niratures  are  the  reconmended  operational  Units 
for  each  Insulation  class  to  Insure  reasonable  lives  in  the  range  of  20,000  to 
40,000  hours.  In  order  to  accelerate  testing,  test  teiaperatuTos  are  selected 
above  the  mlnlasin  recomnMmded  rated  teaperetures.  This  tends  to  force  the 
failure  node  to  one  of  ehasdcal  degtadatlm  of  die  Inaulatlon.  At  tempera- 
tures below  the  recoMamndad  limits,  odter  envlronaantal  and  aeebenlcal  con- 
siderations can  Influence  life.  This  Inolodaa  factors  such  as  dirt,  oil. 
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ch«mlcalt,  oxygen  and  moisture  as  well  as  thermal  stresses  resulting  from 
thermal  cycling.  In  addition  to  environmental  factors,  the  actual  quality 
in  manufacture  of  production  motors,  such  as  nicks  in  the  insulation  durlrtg 
winding,  or  Improperly  seated  wedges  in  slots  can  lead  to  early  failures. 
Quality  in  the  very  small  fractional  horsepower  motors  become  even  more 
Important  because  of  physical  size.  All  of  these  factors  tend  to  Induce  oc- 
casional failures  at  low  temperatures  which  would  not  be  predicted  by  the 
Arrhenius  equation.  It  is  anticipated  therefore,  that  data  from  the  failure 
bank  will  contain  some  early  failures  and  will  not  be  completely  compatible 
with  literature  predictions. 

The  following  equations  for  life  may  be  derived  for  the  average  slope 
and  Intercept  from  the  previous  literature  tabulation  of  winding  life: 


Class 

A (105°C) 

Log 

life 

- 4400/T  - 7.49 

(30) 

Class 

B (130°C) 

Log 

life 

- 5400/T  - 8.78 

(31) 

Class 

F (150°C) 

Log 

life 

- 7000/T  - i2.28 

(32) 

Class 

h (180°C) 

Log 

life 

- 7800/T  - 11.88 

(33) 

These  equations  indicate  motor  insulation  life  including  humidity  cycles. 

Lives  of  3 to  4 times  these  values  can  be  anticipated  when  applied  to  actual 
motors.  As  previously  indicated,  ether  winding  failures  may  be  anticipated 
in  fltotors  due  to  quality  factors  and  environmental  influences.  This  Is 
especially  true  in  the  case  of  fractional  horsepower  motors  since  the  small 
wire  sizes  can  readily  be  damaged  or  over  stressed  during  winding.  To  verify 
this  condition,  the  failure  data  bank  was  reviewed  for  winding  failure  for 
motors  with  stator  laad.natt«in  dlamaterc  from  1.5  to  3.0  inchaa.  The  data  was 
revlswed  consldarlng  the  nunber  of  winding  failures  versus  the  totcl  number 
of  failures.  (The  remaining  fallurea  wars  all  baarlng  failures.)  A similar 
review  was  conducted  on  motors  with  stator  laminations  of  ons  Inch  diameter. 
Table  14  lusmarlzes  these  results. 

For  both  motor  slsee,  it  Is  noted  that,  tha  fallura  percentage  of  wind- 
ings below  e winding  tenpereture  of  160\  ranges  from  3.5  to  11.1  percent.  The 
average  failure  percentage  of  windings  In  this  range  Is  approximately  8X. 


Ttils  is  much  higher  than  predictions  from  the  literature  would  Imply  for  this 
temperature  range  due  to  degradation  of  Insulation. 


lABLE  14 

MOTOR  VTINDING  FAILURES 

(l.S"  - 3.0"  Stator  O.D.)  (Class  F & H Insulation) 


Winding  Temperature 
(By  resistance) 

Total 

Tested 

Winding 

Failures 

Percent 

Windins  Failures 

75 

31 

2 

6.5 

114 

172 

6 

3.5 

130 

137 

15 

10.3 

160 

208 

19 

9.13 

185 

163 

3 

1.84 

(1.0" 

Stator  O.D.) 

(Class  H Insulation) 

75 

27 

3 

11.1 

100 

75 

5 

6.7 

122 

98 

8 

8.2 

175 

392 

69 

17.6 

190 

10 

5 

50.0 

The  aacoiid  point  noted  from  these  tables  is  that  the  smaller  motors 
have  higher  failure  percentages  of  windings  at  the  higher  temperatures,  than 
the  larger  motors.  In  fact,  windings  become  a small  percentage  of  failures 
in  the  laraar  motors  at  high  temperature.  At  these  higher  temperatures, 
bearing  life  drops  significantly;  therefore,  it  ia  expected  that  winding 
failures  would  not  be  too  Importart.  However,  in  the  small  motors,  the 
hi^er  temperatures  eccelerate  winding  failure.  The  larger  motora  are  wound 
with  magnet  wire  of  0.005  inches  and  larger  while  the  smell  motors  use  wire 
of  0.0035  Inches  and  smaller.  Difficulty  In  maintaining  uniform  atrees  dur- 
ing  winding,  the  compactnaas  of  winding  and  the  small  turn  radius  nacassaty 
for  the  small  wire  era  conaidarad  primary  causes  for  the  premature  failure  at 
high  teaperetura. 

In  predicting  motor  life,  it  appears  that  coaaidaratlon  aust  be  given 
to  motor  sisa  as  wall  as  to  the  random  failures  encountered  in  the  low  tempera* 


tur«  range.  IXirlng  analysis  of  tha  failure  bank  data  these  two  areas  were 
considered. 

3.3.2  Regression  Analysis  of  Winding  Failures 

Table  14  outlines  the  results  obtained  by  reviewing  the  number  of  wind- 
ing failures  in  overall  motor  failures.  It  was  pointed  out  that  winding  fail- 
ures represented  between  3.S  and  10  percent  of  the  total  failures  at  winding 
teiiparatures  below  160°C.  If  we  utilize  literature  predictions  for  winding 
life,  the  results  lead  to  lives  orders  of  magnitude  greater  than  bearing  life. 
This  Implies  that  winding  failures  would  not  be  encountered  in  motors  using 
graase-pack  ball  bearings.  Table  14,  however,  does  not  verify  this  condition. 
Therefore,  it  is  anticipated  that  these  failures  are  not  due  entirely  to  In- 
sulation aging  and  may  be  more  random  In  nature. 

Table  4 Is  the  sussnary  of  winding  failures  from  the  failure  data  bank. 
Inspection  of  the  data  Indicates  that  the  significant  portion  of  failures  . 
occur  at  andilent  ten^iera cures  of  12S°C  and  above.  This  also  follows  the 
trend  of  data  presented  in  Table  14  for  small  diameter  (1.0  Inch)  motors. 
Although  not  indicated  In  Table  4,  all  failures  at  125^0  and  above,  are  Indeed 
small  dlaawter  motors.  Regression  analysis  of  life  versus  temperature  there- 
fore will  be  dondnaced  by  the  lower  lives  at  higher  temperature.  A linear  re- 
gression was  conducted  of  the  data  of  Table  4 entering  the  following  veriables 
temperature;  number  of  polos;  voltaga;  spaed;  winding  size;  insulation  type, 
and  smtor  Cypa.  ' The  analysis  indlcati^  that*  t«|nc>Qrature  and  motor  quality 
rapresented  tha  significant  varlablas  influencing  winding  life.  The  results 
of  the  analysis  indicated  winding  life  could  be  expressed  es: 

Log  life  -^’^.3.52  In 
w 

idiere  q « quality  code 

commsrclel  quality  ■ -.13 

4 

Hll-spac.  quality  * -f.lS  . 
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In  comparing  these  results  with  Equations  ('32)  and  (33)  for  Class  F and 
H Insulation,  there  Is  very  poor  correlctlon.  It  would  appear  that  the  lower 
lives  of  the  small  motors  dominate  the  prediction  model  of  Equation  (34). 

To  determine  the  influence  of  the  high  ambient  temperature  failures, 
winding  life  of  units  operating  at  an  ambient  temperature  below  125°C  was 
plotted  as  a function  of  temperature  on  Figure  6.  Although  the  regression 
line  is  drawn  from  only  five  points  It  does  indicate  a temperature  slope  ap- 
proaching that  of  bearing  life,  l.e.; 


Log  life 

2395 

(winding)  ■ 2.064 

(35) 

^amb 

Log  life 

(bearings)  - - 2.534 

(19) 

stab 

If  the  lo%rer  temperature  winding  failures  are  random  in  nature,  it 
%rould  be  anticipated  that  the  Welbull  slope  would  be  1.  Therefore,  the  mean 
and  standard  deviation  of  the  Welbull  slope  for  all  winding  fallursa  as  well 
as  low  temperature  winding  failures  was  conducted  as  tabulated  in  Table  15. 


TABLE  15 

Winding  Failure  Welbull  Slope 


Mean  Values 

Mean  Std,  Dey, 

All  Winding  Failures  1.961  1.195 

At  or  Above  12S^C  Failures  2.329  1.250 

Less  than  125^C  Fallurekf  1.085  .716 

(Bearing  Failures)  (2.878)  (1.997) 


From  the  Welbull  slopes  of  Table  15,  winding  failures  below  125‘’c 
ambient  have  a mean  slope  approaching  1 indicating  random  type  failures  with 
« constant  fellure  rata.  Failures  of  the  l.O  Inch  diameter  antors  at  125°C 
ambient  have  a mean  Welbull  slope  of  2.329  Indicative  of  veer  out  type  fail- 
ures. The  Influanca  of  life  on  taaq>aretura  of  data  In  the  data  bank  therefore. 
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appears  to  be  Influenced  by  two  different  failure  modes.  Equation  (34)  there ^ 
fore,  does  not  appear  to  be  a reliable  model  for  predicting  temperature-life 
characteristics.  The  wear  out  Equation  (35)  for  motors  larger  than  1 inch 
diameter  appears  to  be  a more  realistic  model  for  winding  Influences.  For 
these  motors,  the  percentage  of  winding  failures  versus  bearing  failures 
appear  to  be  almost  constant  over  the  temperature  range.  If  we  make  the  as- 
sumption that  winding  failures  are  a constant  percentage  of  total  failures  over 
the  temperature  range.  Equation  (35)  can  be  modified  to  accommodate  this  re- 
lationship with  the  slope  of  log  life  versus  temperature  the  same  for  both 
bearings  and  windings. 


Log  life  winding  ■ 

Characteristic  life  ■ 
Weibull  Slope  • 

Log  life  bearings  ■ 

Characteristic  life  • 
Weibull  slope  ■ 


m 

22746  hours 
1.0 
2357 


373 


- 2.534 


6096  hours 

2,88 


4.357 


3.785 


Figure  7 illustrates  a Weibull  plot  of  bearing  life  of  the  equation  above 
and  winding  life  at  a slope  of  2.88  and  1.0,  respectively.  The  cuwes  cress  at 
the  11  percentile  point  Indicating  that  approximately  11  percent  of  the  failures 
will  be  electrical  as  a competing  risk  failure  mcdel.  This  is  slightly  higher 
than  the  average  failure  percent  of  8.0  from  Table  15.  The  plot  docs  illustrate 
the  method  to  be  used  for  pr<odictlng  winding  life.  That  la,  the  average  8 per- 
cent wludlng  failures  Is  used  from  Table  15  as  typical  of  the  nunber  of  winding 
failures  of  motors  versus  bearings.  A Vslbull  slope  of  1 Is  for  winding 
failures.  The  winding  life*,  and  bearing  life  at  the  8 percentile  point  are 
equal;  i.e. 

io,  Uf.^  + |l»»  in  • Log  Uf.j  + Ug 

'"g  1“*,  ■ ^ t:?08  I • 1 
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in  (FAILURE -AGE) 


i,og  life  - Log  life,  + .704 
o w B 


Log  Ufe^  - T 


2357 


- 2.534  + .704 


anb 


2357 

Log  life  » - 1.83 

amb 


(36) 


IrA  a similar  manner,  the  life  of  the  winding  versus  winding  temperature 
was  develop.'^d  as 


2342 

Log  life  - - 1.206  + q 

w 1 


(37'» 


V 

whete:  q is  defined  from  Equation  (34) 

These  two  equations  provide  the  analytical  model  for  random  failures 
encountered  in  windings. 


o 

Wiriding  life  at  ambient  temperatures  of  125  C for  motors  with  stator  dla* 
meters  of  1.0  inch  and  smaller  have  reduced  life  from  Equation  (36)  or  (37). 

A regression  analysis  was  conducted  on  only  small  motors  at  teiqperatures  above 
ambient  with  the  following  results: 

Log  life^  - - 6.328  (38) 

w 

Log  life  - - 8.2  (39) 

w 


Two  aquations  wera  developed  one  using  ambient  temperature  the  other 
winding  temperature.  It  wes  stipulate  in  the  sections  discuesing  bearing  life 
thet  it  Is  often  difficult  to  find  the  bearing  tenpereture  without  meesurements 
or  Informfttlon  furnished  by  the  menufacturer.  Tha  sana  problam  exlata  with  pre- 
dicting winding  temperature.  Table  8 provldea  e aucanevy  of  typical  Boteor  tampara- 


ture  rise  data.  Because  of  the  broad  range  of  temperature  rise  that  could  be 
encountered,  a mean  rise  of  36°C  (see  Figure  2)  Is  reconnended  when  better  In- 
formation  Is  unavailable.  To  provide  further  slmpllcatlon.  It  Is  suggested 
^hat  ambient  temperature  be  used  for  the  winding  failure  model  when  actual 
winding  temperature  rise  is  uncertain. 


To  retain  the  simplicity  of  the  winding  failure  model  using  ambient 
temperature,  it  is  suggested  that  the  Influence  of  quality  and  motor  size  be 
disregarded  for  this  model.  The  more  complex  model  utilizing  winding  tempera- 
ture which  necessitates  more  detailed  information  regarding  motor  characteris- 
tics would  include  quality  and  motor  size.  Equations  (36)  through  (39)  would 
reduce  to: 


a)  Winding  life  using  ambient  temperature: 


Log  life  + Log 
w 


1 1 . 2357  , o, 

l-F^t)  T " 


(40) 


amb 


b)  Winding  life  using  winding  temperature: 


1 23&2 

Log  life^  + log  In  ^ - 1-206  + q^ 

w 


(41) 


where:  <1^  " '*1^  comnercial  quality 

« +.13  Mil.  spec,  quality 

c)  Winding  life  using  winding  temperature  and  motors  with 
■tator  O.D.  2.54  cm  or  lees: 


where  ■■  characteristic  life  (hra''  (42) 


Bquetlons  (40^  and  (41'  present  the  failures  at  any  percentile  (F(t)^ 
while  Enuatlon  (42)  Is  the  predicted  life  fpr  the  characteristic  life  only. 

S2 


The  multiplier  for  life  at  other  percentiles  is: 


At 


At 


winding  temperatures  of  160  C and  greater: 

1 \ 2.329 


life  ■ a 


In 


l-F(t) 


winding  temperatures  of  159  C or  lower: 


life  ■ a 


In 


l-F(t) 


(43) 


(44) 


It  has  been  suggested  that  ambient  temperature  neglecting  quality  code 
(Equation  40)  be  used  for  predicting  winding  life  in  the  absence  of  detailed 
knowledge  regarding  the  motor.  Neglecting  the  quality  code  Introduces  a 35 
percent  error  in  determining  winding  life.  This  Is  equivalent  to  an  error  of 
8^C  In  estimating  temperature.  Winding  temperature  rise  using  Equation  (40) 
Is  an  average  rise  of  36°C  as  shown  In  Figure  2.  In  many  applications » pre*> 
dieting  and>lent  temperature  duty  cycle  Introduces  errors  In  the  range  of  lO^C 
or  greater.  The  quality  code  therefore.  Is  less  significant  In  the  Introduc* 
tlon  of  errors  than  Is  the  temperature  Influence  end  the  use  of  a nominal 
qixallty  value  Is  proposed  to  maintain  the  ease  In  use  of  Equation  (40). 

If  It  Is  known  that  the  motor  runs  considerably  hotter  or  cooler  than  36*^C, 
Equation  (41)  should  be  applied. 


Neglecting  small  motors  of  1 inch  (2.54  cm)  stator  lamination  diameter 
and  smaller  Introduced  error  only  when  ainblent  temperatures  exceed  120*^0.  The 
motor  size  and  ambient  temperature  represents  a very  small  portion  of  the  total 
motor  population.  Therefore,  neglecting  small  motors  at  high  temperature  Is 
felt  warranted  for  e generalised  FHP  winding  failure  prediction  model. 


It  should  be  pointed  out  that  in  the  development  of  failure  prediction 
models  for  both  bearings  and  windings,  the  slnpllfled  equations  are  developed 
for  use  In  MIL-HDBK-217B.  The  approach  Is  to  provide  an  ultimate  FHP 
motor  failure  model  chat  requires  the  minimal  of  Infoimatlon  from  the 
user  to  estimate  motor  life.  It  is  felt  that  the  siapllfled  equations  can 
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perform  that  function  if  the  user  is  aware  that  a more  comprehensive  approach 
is  available  for  unusual  applications  of  temperature,  speed,  load  or  environ- 
ment. 


3.3.3  Influenc'i^  - £ AUitude  on  Bearing  and  Winding  Temperature 


A number  cf  individual  tests  in  the  failure  data  bank  were  run  at  pres- 
sures simulating  altitudes  between  50,000  and  85,000  feet.  Bearing  tempera- 
ture rise  at  these  altitudes  %iere  assumed  to  be  40°C.  Since  each  test  was  a 
single  test  point  and  several  tests  were  run  at  each  temperature,  the  average 
life  at  each  ambient  temperature  plus  AO'^C  was  plotted  on  Figure  8.  Also 
plotted  was  the  predicted  life  from: 

/ 


Log  life  ^ 


log  In 


1 


^l-F(t)  j 


2342 

Tb 


1.91 


(10> 


T 


for  both  the  characteristic  life  (F(t)  - 0.632)  and  the  B-50  or  fifty  percent 
failure  (F(t)  ■ 0.5).  The  B-50  Ufa  Is  more  suitable  a comparison  with 
avaraga  life.  It  will  be  noted  from  Figure  8 that  altitude  operation  of 
50,000  feat  introduces  an  undarestlmata  of  winding  temparature  of  appronl- 
BMtaly  95^C.  At  a 125^C  assusiad  bearing  temperature,  the  reduction  In  life 
can  ba  saan  to  ba  10  to  1. 


Unfortunately,  data  is  net  available  at  altitudes  below  50,000  feet 
to  provide  an  Indication  of  the  Influence  of  life  versus  altltutde.  Since 
the  Ufa  influence  Is  so  significant,  a derating  factor  should  be  applied  to 
life  with  operation  at  altitude. 

Vainott  (21)  does  provide  some  indication  of  the  Influence  of  altitude  on 
uotor  temperature.  Figure  9 taaen  from  Reference  21  Indicates  the  temparature 
rise  eharpctprlstics  of  intamelly  fan  cooled  and  non-fan  cooled  totally  en- 
olosad  motors  versus  altltutde.  At  altitude  lavalo  of  50,000  feat,  the  tampar- 
atura  rise  of  fan  cooled  aotera  are  approxlwctaly  5 ttmos  the  saa  level  temper'' 
atura  rtaa  while  In  totally  ancloaed  motors  Lha  rite  is  2 to  7.  that  at  aaa  lev- 
el. Thla  anplalna  the  higher  iMtarlng  tamporatura  rises  aasuclatad  with  Figure 
8.  Thaaa  nushera  howavar,  ara  aosarahat  mis leading  tinea  the  ambient  tampara- 

S4 


Flgur«  9 

T««p«riitur«  Charaotariitica 
Internally  Fan  Cooled  and  Non>Fan  Cooled 
Totally  Bncloaed  Motora 
Varaue  Altitude 
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ture  also  decreases  with  altitude.  At  50,000  feet  for  standard  conditions, 
the  ambient  temperature  Is  >55°C.  For  a fan-cooled  motor  having  a 50^0  rise 
at  sea  level,  winding  temperatures  can  reach  200°C.  At  sea  level,  the  winding 
temperature  would  be  65^C.  The  Influence  of  50,000  feet  altitude  therefore 
would  be  3 to  1 over  sea  level  temperature.  In  a totolly  enclosed  motor,  at 
50,000  feet,  the  winding  temperature  would  reach  50*^0  compared  to  65°C  tai^>ar- 
ature  at  sea  level.  This  simply  points  out  that  for  a fan'-tooled  motor  with 
a sea  level  winding  temperature  rise,  altitude  will  not  result  In  higher  wind- 
ing temperatures  below  35,000  feet.  Above  this  altitude  rapid  reduction  In  life 
can  be  anticipated.  For  totally  enclosed  motors,  altitudes  to  60,000  feet  will 
not  influence  life  as  compared  to  sea  level.  In  each  case  it  Is  assumed  the  am- 
bient temperature  follows  the  standard  temperature-altitude  characteristics. 


From  these  results,  the  influence  of  altitude  on  life  for  a fan-cooled 
motor  may  be  tabulated  baaed  on  a 50  parcent  reduction  in  life  for  every  lO^’c 
increase  in  tenperature.  Table  16  is  a tabulation  of  life  multipliers  versus 
altitude  based  on  Army  maxinhim  condition  temperatures.  « 


TABLE  16  (Reference  21) 


Influence  of  Altitude  on  Motor  Life  for  Pan-Cooled  Motors 


Altitude 

(thousands  of  feat) 


Life  (kiltipliers  for 
Sea  Level  Motor  Temperature  Rise 

2Q**C  _ 30**C  _ _ A0°C 50°C  60°C 


S.L. 

1.0 

l.O 

1.0 

l.O 

l.O 

25 

1.0 

1.0 

1.0 

l.O 

1.0 

30 

1.0 

1.0 

l.O 

1.0 

0.5 

40 

1.0 

1.0 

0.25 

0.0078 

0.0625 

50 

1.0 

0.125 

0.0078 

0.0005 

60 

0.06 

0.0005 

It  is  recommended  thet  the  Mtor  life  be  ooeputed  tot  sea  level  condi- 
tions and  22^C  using  equatloas  of  Table  17.  Lii^  Is  then  derated  per  Table  16. 


Alternately,  the  bearing  tesperatiure  and  viatfiitg  temperature  «an  be  eatUpted 
fron  lable  8 and  appllad  to  tho  oguationa  of  Tabla  17. 
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3.4  FHP  Motor  Follurt  Modal 


la  SactioTM  3.2  and  3.3  aquatlona  iiara  davalopad  to  pradlet  a baarlng 
failura  Modal  and  a winding  fallura  nodal  raapactlvaly.  Ihaia  two  fallura 
Bodat  doodnata  tha  ovarall  Ufa  of  FHP  notora.  Slnca  aithar  fallura  may  occur 
during  tha  Ufa  of  a notor  thay  ara  conaldarad  aa  conpatlng  failure  modes.  In 
Pafaranca  1 a conpatlng  fallura  nodal  was  davalopad  basad  on  tha  following  dls* 
trlbutlon: 


F(t)  • 1 - a 


(45) 


whara:  F(t) 
t 


fallura  percantlla 

tlna  to  fallura  (hrs) 

baarlng  charactarlstlc  Ufa  (hrs^ 

Walbull  slopa  for  baarlng  falluraa 

winding  characteristic  Ufa  (hrs) 

Walbull  slopa  for  winding  falluraa 


By  caking  tha  natural  Icgarlthn,  Equation  (45)  nay  be  rewritten: 


ln(l-F(t)) 


a 

w 


P. 


To  solve  chis  aquation  for  notor  Ufa  (t)  knowledga  of  tha  Walbull  slopa  and 
charaatnAatlc  life  nust  ba  obtained.  Equations  (26)  and  (27)  describe  tha 
fallura  nodal  for  baailnga  providing  the  value  of  o.  Equation  (26)  requires 
considerable  knowledge  regarding  the  notor  and  bearing  while  Equation  (27) 
aaenssltataa  knowledge  of  the  aableut  tenperature  only  ae  the  sacrifice  of 
oppronlnetely  11  percent  In  the  accuracy  of  the  prediction.  Figure  5 is 
a distrlbntlon  of  Walbull  slopes  indicating  a naan  slope  of  of  2.878.  It 
wilt  be  reeegnlaed  linMdlately  that  solution  of  Equation  (^)  will  require 
Itevatlssi  teehnlques  slnee  la  raised  to  the  2.878  power.  Changing  the 
aquation  to  a cubic,  whleh  Is  within  one  standard  devlatlen*  slnpllflas  tha 
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solution  of  the  equetion.  In  the  general  form  therefore,  a value  of  « 3 le 
suggested  while  Pg  ■ 2.678  is  suggested  for  the  more  complex  solution. 

Windl.xg  characteristic  life  can  be  computed  from  Equations  (40), 

(41),  or  (42).  Equation  (40)  Is  the  simplified  general  winding  life  model 
using  ambient  temperature.  Equations  (41)  or  (42)  are  the  complex  form  In- 
cluding motor  slae,  quality  and  winding  temperature.  The  Welbull  slope  for 
winding  £allur«,s  as  outlined  In  Table  15  Is  p^  « 1.0  for  all  motors  except  those 
I Inch  (2.54  cm)  outer  diameter  of  the  stator,  operating  at  winding  temperatures 
of  IPO^C  or  greater.  For  this  special  case  P^  ■ 2.329. 

Table  17  summarises  the  overall  FHP  motor  failure  models  from  Equation 

(45), 


TABLE  17 

FHP  Motor  Life  Prediction 


A.  General  Case  (simplified  Model) 


A-l 

where: 

In  (l-F(t)) 

1. 

• ■ ■ 

1 

*■  H 

wJ 

1 

A-2 

“a  • 

10 

_ 2.534 

^amb 

+ 

10 

'-fSS.  +20 

amb 

s « 

A-3  bog  ^ - 1-83 

Amb 


(eemtinued) 


59 


TABLE  17  (continu«d  from  prtvioas  p«g«) 
Complmx  Casa  (coinplota\ modal) 


(Equation  B-1) ; Inil-F(t) 


2.878 


(Equation  B-2) : ln|l-F(t)|  ■ - 


2.878 


whara: 


for  motors  greater 
than  1 In  (2.54  cm) 
stator  outer  diameter 
and  all  motors  with 
winding  temperature 
less  than  159°C. 


2.329 


for  motors  equal  or 
less  than  1 in.  (2.54 
cm)  stator  outer  dia- 
meter at  winding  ^ 
temperature  of  160  C 
or  greater 


(Equation  B-3) : 1 


+ . - 4.32  D»  X lo'*  - It,  - .oomj  Ip  1 ' + 


+ 1,. 7 

10  »■  ® J + 300 


(Equation  B-4):  Log  " 1»206  + q^  (for  use  with  Equation  B-l) 


(Equation  B-5) : 


riM 

T. 


1.206  +^  + r450< 

10  L'*' 


4500  - 6. 328-0 J 


(for  use  with  Equation  E-2) 


It  is  obvious  from  Inspection  T&ble  17  that  the  simplified  model 
offers  considerable  advantage  in  the  ease  of  use.  Table  6 provided  an  estlmata 
of  the  accuracy  of  prediction  using  the  complex  model.  It  is  noted  that  75 


m: 


parcant  of  the  failures  would  be  predicted  within  e life  model.  Althou^ 
the  prediction  accuracy  is  improved  11  percent,  the  ability  to  estimate  all 
constants  in  the  complex  model  decreases  the  probable  prediction  accuracy. 

Ihe  simplified  oiodel  therefore,  was  selected  as  the  general  prediction  model 
for  use  in  the  reliability  handbook  with  stipulations  regarding  its  suitability 
for  extremes  of  temperature,  speed,  load  and  motor  sice.  For  those  type 
applications  the  more  cosiplex  model  is  recoomended. 

3.4.1  Failure  Rates  of  FHP  Motors 


The  method  isost  often  used  to  analyte  failure  populations  is  thru 
detersdnatlon  cf  the  failure  rate  of  the  population.  The  failure  rate  also 
termed  the  hasard  rate  over  the  interval  e^«t:ft^  ** 

ratio  of  the  nuieber  of  failures  occurring  in  the  time  interval  to  the  number 
of  survivors  at  the  beginning  of  the  tisw  interval  divided  by  the  length  of 
the  tism  interval;  i.e.: 


.(t.)  - 100  j 


The  hasard  rate  is  a measure  of  the  instantaneous  rate  of  failure  and  Is 
useful  in  describing  the  type  of  failure.  Increasing  failure  retes  signify  a 
wear  out  process;  constant  failure  rates  indicate  a random  failure  mode  and 
decreasing  failure  rates  are  associated  with  infant  mortality  failures  gener- 
ally due  to  design  weakness  or  assembly  defects. 


For  Increasing  failure  retes  such  as  encountered  with  bearings,  the 
instantaneous  failure  rate  at  any  tine  (t)  ia  not  a good  estimate  of  the 
numbsi*  of  failures  that  have  occurred  to  tiae  (t).  The  everege  failure  rate 
H (c)  of  the  aotor  at  the  SO  purceatile  failure  is  the  suggested  failure  rate 
to  be  coabined  with  tha  other  constant  fetlure  tste  coaponents.  Tha  avaraga 


failura  rata  can  ba  dafinad  as: 


• D 


h(t)  dt 
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If  the  average  failure  rate  to  any  time  (t^)  from  time  0 Is  utilized,  Equation 
(46)  reduces  to 


H(t) 


avg 


h(t)  dt 


However,  the  integral  of  h(t)  is  the  cumuJative  hazard  function  H(t).  There- 


fore- 


For  the  con^etlng  failure  model,  the  cumulative  hazard  function  was  defined 
as : 


B 3 

H(t)  - (— ) (“) 

B w 


(47) 


The  average  failure  rate,  therefore,  is 


H(t) 


. KiJCl 


avg 


H(t) 


.1^ 

evg  P 

a 

B 


P -I 

V 


(48) 


a 


Table  17  defines  the  values  of  and  P„*  The  time  (t)  is  the  »• 

g if  B V » ; 

nocor  life  as  calculated  from  the  equations  of  Table  17.  If  we  consider  only  | 

I 

the  alible  case  of  Table  17,  Equation  (48)  can  be  rewritten  as: 
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H(t) 


avg 


^ + ^ X 

“b  . 


10^  failures/ 10^  hours 


where 


‘B 


as  — 

|”57  .2.534) 

-4500  U 20 

10 

[ amb  j 

10 

/ + 300 

<48a) 


In  l-F(t) 


a.  - F - 

anb 

and  t is  calculated  from: 
" 3 

a + a ! 

B w J 


(49) 


The  mean  time  to  failure,  as  previously  discussed,  is  estimated  at  the 
50  percentile  failure  condition  where  F(t)  « 0.5.  The  natural  log  of  0.5  is 
-.69315.  Equation  (49)  may  be  rewritten  as: 


.69315 


This  is  a cubic  in  t and  may  be  solved  for  t. 


(50) 


t - 


34657  +V.  12011  + . 

L 


03704 


\l  I 

V. 34657  -V.  12011  + ,03704 


3.4.2  Use  of  FHP  Motor  Failure  Models 

In  order  to  assist  the  reader  in  the  use  of  FHP  motor  failure  models 
sad  prediction  techniques*  so  exuple  has  been  selected  that  utilises  the 
developed  equations.  loth  the  simplified  model  and  the  complex  model  are 
applied  to  illustrate  predictions  using  both  procedures. 

63 


It  It  detirable  to  obtain  motor  life  and  failure  rate  of  a 1/4  HP 
60  He  motor  that  operates  50  percent  of  the  time  at  an  18°C  aad>lent  and  the 
remainder  at  lOO^C  aiid>lent.  Ihe  motor  Is  assumed  to  be  of  connmrclal  quality, 
standard  end  bell  construction,  a capacitor  start  Induction  motor  vlth  silicone 
lubricated  beatings  that  are  loaded  with  a ten  pound  load  per  bearing. 


3.4.2. 1 Simplified  Model  ^ 

I 

Id  this  model.  It  Is  necessary  to  use  only  the  ambient  temperature  for  | 

prediction  of  motor  life  and  failure  rate.  In  the  case  cited,  operation  occurs  | 

at  two  different  temperatures  and  bearing  and  winding  life  must  be  computed  | 

for  both.  From  Table  17  Equation  A-2  end  A-3  are  used. 


a)  Bearing  Life  (Eq.  A«2) 


i 

I 


lim  . 2.534 

10 


^ -20 

10.  + 300 


1 - .000164  + .000000011 


* 6095  hri 


Overall  bearing  life 


Equation  (21)  Is  used  to  obtain  life  under  thermal  cycling  conditions. 


^1  , 111 


where  N 


• no.  cf  cycles,  h^  • hrs.  at  T^,  and  Og^  ■ life  at  T^. 


By  substituting  -jpVlT  **  »idiere  Og  ■ cycling  characteristic  life. 


Chen: 


^ V 7*2  / “b  Y 1 

* **2  \ j ^2  \ % / 


Since  proportion  of  time  at  each  temperature  ■ .5(or  SOX),  we  have: 


“b  + •*  “c 


.5  og  . -So, 

TTWr  * ■JBB  ' 


10ZZ3  hr. 


b)  Vlndins  Lit*  (Bq.  A-1) 
,r  iflo  Aablfflt 

Log  a 


ilSL 


UB3 


1 0 \ ■ - . 

a - 1860614  hrt 

*1 


|00  C Ambl«i£ 

“ w “ ^ ‘ 


a - 30834  hr« 


Ovorall  Llfo  (Equation  21) 


.5  a 


.5 


1860614 


30834 


- I 


a ■ 60663  hrs 


c)  Hotor  Ufa  (Equation  A-1,  Equation  49) 


. . I' 


• ln(l-r(t)) 


at  F(t)  - 0.^ 


1-3  \ ^ ‘ 3 , — ^ “ 

t - 10223|Y  .693394  + -\^.000255 

: j j j -■  V ' ’ ■V.'Jj'"'  ‘ -**'■ 


t • 10223  I .0651 


[■■ 


- .O6337I 


i 


(A-1) 


r f' 


t • 8400  lira 


d)  Motor  falluro  rota  (Equation  48A) 


t I 10223-*  ® 


,]xl0‘ 


- 82.3  falluras/io''  hr« 


3. 4.2. 2 Coaplax  Modal 

In  this  Bodal  It  la  naea^sary  to  obtain  winding  toaparatura,  baarlng 
taaparaturo  and  baarlng  cliaractarlatlca.  Vo  obtain  tha  oparatlng  taoparatum 
rafar  to  Tabla  8.  For  tha  not or  In  quaatlon,  tha  high  slip  eatagory  la  uaad 

with  taa^MraCura  riaaa  aaVtctad  at  low  allp;  l.a.  motor  tanperatura  riaa 
30^C,  baarlng  tawparatura  rlae  15*’c.  For  aora  accurata  raaulta  tha  nanufaeturar 
should  ba  conaultad. 

a)  Baarlng  Ufa  (Tabla  17  Equation  B-3) 

It  la  nacaaaary  to  aatabllah  quality*  graaaa  typa*  DN  and 
load  aa  followa: 


quality  (q):  Saa  Equation  <17);  coaawrclal  q ■ -.27 
graaaa  (X  ) : Saa  Sactlon  3.2. 2.6  graaaa  5;  • 1.81 

m : Saa  Tabla  9 

.23HP  ■ 186  watts;  D ■ 13  m 

Equation  (20) : M ■■  1^  » 60  ■ 3600  RFM 

2 

DM  « 46800 


SP  for  13 
HI  bMring; 
(Specific 
load  at 
33-1/3 
r?«) 


glvan  as  10  Iba 

Saa  brg.  uifgr.  catalog  for  0.5  In.  brg. 
SP  - 505  (R-8  brg) 


18  + 15  273  • 306°K 

and 

100  + 15  + 273  - 388*^C 


8ra.  Ufa  at  18°C 

Froa  Equation  (18) 

Log  Ufa.  - + q - 4.32  DN  x lO"® 

p 

- - .27  - 4.32(.0468)  - 1.81  - 3.6j^^j^*^ 

- .00000435  + .000070 

- 13414  hrs 

Bra.  Ufa  at  10Q**C 


Tw 

- K - .001  Mj  — 


Bcpiatlon  (18)  partially  solvad  abova;  l.a. 

Log  Ufa  • - 2.293 

Log  Ufa  ■ 3.743 
Equation  (1-3)  as  abova 


X 


+ ■■  - -11^ 


*2'  10^*^*^  lOj^  388 


^P  + lT? 


* 300 


Bj  - 5539  hrs 
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O^fTfU  lfU,t  Cgwa£teB-aU 


.5  a 


B 


B, 


.5  0- 


B„ 


- 1 


• 5 a 


.5 

13414  5539 


B 


- 1 


Ojj  ■ 7840  hri 


b)  Winding  Life  (Table  1/  Equation  B»4) 

Log  “ - - 1.206  < <! 

^ ^wlnd 


T«  - 18  + 30  + 273  - 321  L 


- 100  + 30  + 273  - 403  K 


From  Equation  (42)  q^  ■ .13 


Life  at  18”C 


Log  - 1.206  - .13 

a - 311906  hrs 
’'l 

igp!g. 

a • 29882  hra 

’'2 

Overall  Wladloa  Ltfli  lBqMtigB  2D 

.5  .9 

9iiw  SfSM  " ^ 

- 57868  hrs 
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c)  Motor  Llfo  (Bn-  Ll*o)  Toblo  17  Equation  (B-l) 


ln(l-F(t)) 


■'h 


2.878  ^ '] 


+ 


^1 


ln(l‘^).!')  -0-693147 

2 878 

-693W  ^ ‘ + 5^ 


Tfela  aquation  vaa  aolvad  by  computar 
t - 6492  hra 


d>  Motor  Fallur#  Bata 
r. 1.878 


r i.a: 


10^  failuraa/lO^  hra 


5492__1»878  g 

7840  2-878  57868  ^ 


X • 107  falluraa/lO  hra 
t 


Tha  raaultc  of  tha  complex  model  indicate  a failure  rate  X ^ of  107 
mhila  the  aiapla  model  X^  la  82. S. 
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4.0  CONCLUSIONS 


A review  of  failures  and  failure  nodes  of  FHP  motors  was  conpleted  and 
a reliability  prediction  model  developed.  The  following  specific  conslusions 
and  recosnmndations  ware  developed  as  a result  of  these  studies: 

1.  Two  failure  modes  dominate  the  life  of  FHP  motors  utilising 

grease  packed  non-relubricated  ball  bearings;  i.e.  bearing 
failures  and  winding  failures.  Bearing  failures  represented 
approximately  80  percent  of  the  total  failures  of  these  type 
motors.  Winding  failures  were  more  prominent  in  the  early  life 
of  the  motor  and  were  attributed  to  shorts  and  grounds  resulting 
from  assembly  quality  rather  than  long  tens  insulation  degrada- 
tion. The  winding  failure  percentage  of  soiall  motors  of 
2.54  cm  lamination  OD  operating  at  aBd>lent  temperatures  of  125*^C 
increased  to  50  percent  which  was  attributed  to  the  small  wire 
size  used  in  these  motors. 

2.  Both  failure  modes  were  analyzed  utilizing  a Welbull  cumulative 

distribution  function.  A Welbull  shape  parameter  in  the  range  of 
2.8  to  3.0  was  used  to  describe  bearing  failuras  indicating  an 
increasing  failure  rate  and  wear  out  condition.  Winding  failures 
were  described  by  a Wbibull  shape  parameter  of  1.0  indicative  of 
a constant  failure  rate  and  random  type  failure. 

3.  Mil ti -variable  linear  regression  techniques  were  utilised  to 

develop  a mathematical  a»del  of  both  failure  modes.  Life  of 
the  bearings  and  windings  was  most  sensitive  to  temperature  of 
operation.  A simplified  model  was  developed  for  both  failure 
sK)dea  as  a function  of  asdiient  teaqiereture  only.  This  model 
proved  sufficiently  accurate  for  general  predictions  of  motor 
life.  Secondary  influences  including  SK>tor  load,  speed,  grease 
type,  winding  temperature  use,  bearing  temperature  rise  and 
quality  were  used  to  develop  a more  oomplaa  model  for  life  pre- 
diction. This  model  was  suggested  when  unusual  operating  con- 


ditlQOA  xMira  anticipated.  Tha  added  cooplaxity  of  this  model  was 
not  daaoMd  naeaaaary  for  most  lightly  loaded  FHP  motor  applica- 
tions. 

4.  Tha  two  failure  modes  uere  raeosribined  as  a competing  riak  failure 

model  for  predicting  overall  motor  life.  Equations  were  also 
developed  for  prediction  of  the  average  failure  rate  vising  the 
50  percentile  failure  point  as  an  estimate  of  the  mtiaa  time  to 
failure. 

5.  Appendix  B was  developed  as  a recoamended  '* failure  prediction 

nodal"  replacesMnt  for  both  fractional  horsepower  motors  and 
electronic  cooling  devices  in  MIL-HDiC-217B. 


5.0  REFSRBMCBS 

1.  Wilson,  D.S.,  Wilkinson,  J.  W. , "Rsllsblllty  Modsl  for  Ministurs  Blow- 
ers,” RADC-TR-7 5-178,  March  1975. (AD  A013735) 

2.  Helson,  W. , ''Hazard  Plotting  Methods  for  Analysis  of  Life  Data  with 
Different  Failure  Modes,”  Journal  of  Quality  Technology,  Vol.  2,  Mo.  3, 
July  1970. 

3.  Poon,  S.  Y.,  "An  Experlnental  Study  of  Grease  in  Elastohydrodynamic 
Lubrication,”  ASME  Paper  No.  72-LUB-B,  1972. 

4.  Scarlett,  N.  A.,  "Uae  of  Grease  in  Rolling  Bearing.!,”  Proc.  Instn. 

Mach.  Engrs.,  1967-68,  Vol.  182. 

5.  "Lubrication  Life  Forecast  and  Teat  Procedure,”  New  Departure  Handbook, 
P.P.-IM-3/64,  Cr^e  GL. 

y 

6.  Bushier,  F.  A.,  Cox,  D.  B.,  etc.,  "Reaearch  Studies  on  Lubricating 
Grease  Compositions  for  Extreme  Environments,”  NLGl  Spokesman,  October 
1964. 

7.  Crisp,  J.  N.,  Wells,  E.  E.,  "Low  Temperature  Performance  of  Greaaea 
in  Railway  Roller  Bearings,”  Journal  of  ASLE,  July  1963. 

8.  Demorest,  K.  E. , McMurtrey,  E.  L.,  "An  Evaluation  of  Grease  Type  Ball 
Bearing  Lubricants  Operating  in  Various  EnvlronMnts,”  Status  Report 
Mo.  2,  George  C.  Marshall  Space  Flight  Center,  1976. 

9.  Christian,  J.  B.,  "Lubricant  Variables  Influencing  Miniature  Bearing 
Torque,"  A.  F.  Materials  Lab  Report  AFML-TR-76-135.(AD  B013912) 

10.  Mathes,  K.  M.,  "Aging  of  Ssiall  Motor  Insulation,”  AIEE  Trans.,  Vol. 

71,  Part  3,  Power  App.  and  Systema,  pp.  2S4i’262,  1952. 

11.  Beebe,  T.  J.,  "Insulation  and  Temperature  Rise,”  Machine  Design, 
pp  9-10,  April  9,  1972. 

12.  Whitman,  L.  C.,  Doigan,  P.  "Calculation  of  Life  Characteristics  of 
Insulation,"  AIEE  trans.,  Vol.  73,  Part  I.,  Cosasunicatlon  and 
Electronica,  pp.  193-198,  1954. 

a 

13.  Barberlch,  L.  J.,  Dakin,  T.  W.,  "Guiding  Principles  in  the  Thermal  • 
Evaluation  of  Electrical  Insulation,"  AIEE  Tech.  Paper  No.  56-248, 

1956. 

14.  Grant,  G.  at  al.,  "Motor  Tasts  Evaluate  Thermal  Enduranca  of  Class 
K Insulation  and  Silicone  Varnish,"  AIEE  Trans.,  Vol.  68,  Part  2, 
pp  1133-1138,  1949. 


73 


15.  Horton,  w.  H. , "A  Statistical  Method  for  Predicting  Insulation  Life 
froa  Sxperiaental  Data."  AIEE  Tech.  Paper  No.  56-177,  1956. 

16.  Hezaan.  C.  J.,  "Motor  Insulation  Life  as  Measured  by  Accelerated  Tests 
and  Dielectric  Fatigue,"  AIEE  Trans.  Vol.  72.  Part  3,  Power  App.  and 
Systssia.  pp  986-991,  1953. 

17.  Leape,  C.  B.  et  al.  "A  Method  of  Evaluating  Insulation  Systems  In 
Motors,"  AIEE  Trans.,  Vol.  72,  Part  3,  Power  App.  and  Systems,  pp  793- 
799,  1953. 

18.  Calr,  T.  J.,  "Insulation  Systems  for  Random-Wound  Motors  Evaluated  by 
Motorette  Tests,"  AIEE  Trans.,  Vol.  73,  Part  3B,  Power  App.  and 
Systesui,  pp.  1702-1707,  1954. 

19.  "Rcflneaents  lA  IFEF  117  'est  Procedure  for  Evaluation  of  Life  Ex- 
pectancy of  Random  Wound  Kotor  Insulation  Systems,"  IEEE  Paper 
No.  68TP636.  1968. 

20.  Brancato,  E.  L.,  Johnson,  L.  M.,  "Thermal  Aging  Probleas  ir  Electrical 
Insulation,  Past,  Present  and  Future,"  IEEE  Paper  Nu.  68C6-EI, 

8th  Electrical  Inaulation  Conference,  December  9-12,  1968. 

21.  Veinott,  C.  G.,  "Effect  of  Altitude  on  Ventilation  and  Rating  of 
Aircraft  Electric  Machines,"  Electrics)  Engineering  Trans.,  February 
1946,  Vol.  65,  pp.  84-90. 

22.  Dowson,  D. , Hlgginson,  G.R. , "Blaato-Hydrodynamic  Lubrication," 

Psrgsarm  Press  Ltd.,  London,  1966.  pg.  36. 


•v 


74 


6.0  LIST  OF  SYMBOLS 


A 

B 

°C 

c. 

c 


D 

D.F. 

DN 

DN^ 

e 

E 

f 

F 

£(t) 

F(t) 

G 

h 

hrs. 

h(t) 

h 

H(t) 

i(t) 

Hz 

o„ 


avg 


Constant 

Constant 

Denotes  bearing  failure 
Temperature,  degrees  centigrade 
Characteristic 

Denotes  combined  bearing  and  winding  failures 

Total  squared  error  coefficient 
Bearing  bore  diameter  in  millimeters 
Degrees  of  freedom 

Bore  diameter  (mi)  times  speed  (rpm) 

Limiting  DN 

Constant 

Constant 

Line  frequency  hs 

Statistical,  F-test  value 

Statistical,  probability  density  function 

Statistical,  ctnaulative  failure  distribution  function 

Dimensionless  sksterial  paramr;ters 

Time  at  teaiperature  (hours) 

Hours 

Statistical,  hazard  function 
Dimensionless  film  thickness 
Statistical,  cumulative  hazard  fiaction 
Statistical,  everagft  failure  rate 

Cycles  per  second 
Degrees  Kelvin 

Knmiber  of  parameters  in  regreasion 
Grease  constant 

Life  (unless  noted  otherwise,  L Is  the  Uelbull  character let ic 
life  when  F(c)  .632) 

Log  to  base  10  (common  logarlthMi) 

Lower  qualiw/  grade 

Log  to  base  a (2.71828)  (natural  logarithms) 


Log 

LQC 

In 
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ORIG 

P 

q 


R 


RPM 

S.D^ 

s,r. 

SQS 


t 

9 

A 


Multiplier 

Mean  sum  of  square*  of  errors 
Nuaber  of  units  or  tharauil  cycles 

Speed  in  revolutions  per  minute,  or  number  of  data  points 

Original 

Number  of  poles 

Quality  code 

Quality  code  for  winding  failures 
Correlation  coefficient  r 

Square  of  the  correlation  coefficient 
Coefficient  of  multiple  correlation 
Revolutions  per  minute 
Standard  deviation 
Specific  load 
Squares 

Sum  of  squares 
Tine  in  hours 

o„ 

Temperature,  K 
Life  coriection  factor 
Dinensionless  speed 
Upper  quality  grade 
Load  in  pounds 
Denotes  winding  failures 
Independent  variable 

Average  of  all  Independent  data  values 
ith  independent  variable  data  value 
Calculated  data  value 
Dependent  variable 

Average  of  all  dependent  variable  data  values 
1^**  dependent  variable  date  value 
Calculated  Deta  value 

llelt!ell  charecterletic  life 
Uelbull  shape  paranater 
residual  error 
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w. 


r 

I 0 

I 

X 

SUBSCRIPTS 

1 aattk  refers  to  aisblent 

avg  refers  to  average 

g refers  to  bearings 

1 refers  to  1th  data  value 

w refers  to  winding 

jj.  refers  to  independent  variables 

V refers  to  dependent  variables 


Lubricant  viscosity 
Standard  deviation 
Percent 


APPERDIX  A 


MULTIPLE  LINEAR  REGRESSION  ANALYSIS 
TECHNIQUES  & EXAMPLES 


A.  IHTRODUCTIOM 

Thia  appandlx  contains  a genaral  daacrlptlon  of  applied  regression  pro« 
cedurea.  Application  of  these  procedures  is  widely  used  today  by  management 
and  acientlfic  personnel.  Successful  use  requires  an  understanding  of  the 
underlying  theory  as  well  as  the  practical  data  problems  which  can  arise  from 
everyday  applications. 

The  following  subjects  are  addressed: 

Regression  Analysis 

Use  of  Regression  Analysis 

Correlation  Coefficients  and  Matrix  Arrays 

Correlation  Coefficients  of  Transformed  Data 

Selecting  an  "Optimum  Set"  of  Explaining  Parameters 

A Large  Number  of  Independent  Variables 

Regression  Tables  Explained 

Examples  of  Data  Reduction 

A.  1 R«3f ression  Analysis 

Regressloa  Analysis  is  a technique  which  can  be  used  to  describe  and 
establish  relationships  between  two  or  more  sets  of  data.  It  is  particularly 
useful  in  cases  where  the  data  is  statistical  in  nature  and  has  scatter  if 
plotted  graphically.  The  relationships  derived  from  this  technique  may  not 
always  provide  a perfect  description  of  associated  data;  nevertheless,  it 
cam  be  a powerful  tool  even  when  approximations  to  data  variatituis  are  fo^md. 
Problems  that  could  only  be  handled  by  guesswork  csn  be  confidently  handled 
once  Imterrelatlonshlps  are  established  with  regression  techniques. 

A, 2 Ose  of  Ragrassloo  Analysis 

Rsgrsssioo  procedures  can  be  used  tot 
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a) 

b) 

c) 


describe 
control,  or 
predict 


In  each  case  one  must  use  a set  of  data  which  Is  numeric  in  form.  In 
some  Instances  the  data  can  be  separated  Into  either  dependent  or  Independent 
groups.  The  dependent  response  variables  are  usually  what  Is  to  be  described, 
controlled,  or  predicted.  In  other  cases,  a knowledge  of  which  variables  are 
"independent"  are  not  known  before  examining  the  data.  Regression  procedures 
can  be  employed  to  establish  the  interdependence  of  such  sets  of  data. 

In  some  fields,  engineering  concepts  or  models  can  be  used  to  help 
select  the  independent  variables  to  be  employed  in  the  regression  relation- 
ship. Once  established,  the  functional  form  of  the  relationship  can  be 
used  to  predict  or  control  the  system  response  from  a given  set  of  independrat 
variable  values. 


llie  power  of  regression  analysis  Is  best  used  whenever  the  response  of 
the  system  to  be  described  is  affected  by  several  operating  parameters  simul- 
taneously. Often  times  system  response  to  a single  independent  operation 
parameter  is  uot  available.  Either  single  variable  control  tests  on  the  sys- 
tems are  much  too  costly  or  cannot  be  Implemented  In  practice.  Given  a col- 
lection of  d^ta,  regression  procedures  can  provide  Insight  as  to  which  varia- 
bles, If  any,  are  most  useful  in  describing  the  system  response. 


A. 3 Correlation  Coefficients 


The  association  between  two  variables  Is  often  examined  by  plotting  In 
graphic  form  one  variable  against  the  other.  Whenever  a large  number  of 
variables  arc  involved  thla  process  la  not  always  practical.  Tha  calculation 
of  a correlmtlom  coafficlamt  Is  one  screenlag  tecaal^pie  idiich  avoids  data 
plottiag  but  cam  still  provide  lnf<»Batlom  on  data  intardapendanca.  The 
"corralatlom  eoafficlamt"  daflnad  by  (R)  balow  Is  one  measure  of  the  dagraa 
of  asaoelatlon  betwesm  two  variolas. 
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The  values  represented  by  S end  Y ere  the  Mens  of  all  the  and  values 
Whose  association  Is  being  tested.  The  o » o are  the  standard  deviations 

^ y 

of  all  X^  and  Y^  values.  A positive  or  negative  (R)  results  whenever  the 
slope  of  the  association  between  variables  Is  positive  or  negative,  respective- 
ly. 


It  can  be  shown  that  the  correlation  coefficient  (R)  Is  always  within 
tha  range: 


-1$  R^  1 

A ssnple  set  of  correlation  coefficients  Is  shown  In  Figure  A-1  with  n dis- 
play of  the  graphic  data  used  to  obtain  the  coefficients.  As  shown,  the 
closer  the  (R)  value  Is  to  the  extrceie  values  of  (+1)  or  (-1)  the  less 
scatter  there  la  in  the  data.  As  (R)  approaches  sero  the  scatter  approaches 
a randan  pattern  for  the  variables  being  tested. 

A. A Correlation  Matrix 

Correlation  coefficients  are  often  calculated  for  more  than  just  me 
pair  of  variables.  In  fact  for  any  nusber  of  (k)  variables  there  can  be 

2 

unique  pairs  of  linear  associations  whidi  night  be  evaluated.  Correlation 
coefficients  are  often  displayed  In  array  or  nstrix  fom  for  convenience.  The 
elenents  of  the  array  are  the  <R)  values  cnlculated  fron  each  pair  of  variables 
whose  association  Is  being  tested.  The  row  and  coluen  nunb«ire  of  the  natrlx 
nerely  represent  the  variables  for  idiich  the  eorre.^ation  coefflelents  have 
been  calculated.  For  exasiple,  see  Figure  A-2.  The  last  colunn  of  that 
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figure  Is  all  (R)  values  for  the  Interdependence  between  teiapereture  and  the 
characteristic  life  of  the  FHP  motor.  The  row  and  column  variable  asalgn^ 
ments  are  as  tabulated  at  the  bottom  of  the  figure.  The  first  rot^column 
entry  (.99)  is  a number  very  cloae  to  (1.00)  and  represents  the  high  correla- 
tion of  the  ambient  temperature  data  with  Itaelf.  All  diagonal  elements  of  a 
correlation  matrix  will  be  very  near  1.00  If  a linear  association  between 
variables  is  being  tested.  In  Figure  A-2,  the  diagonal  alcmentn  are  all  one 
or  near  one.  but  are  shown  as  (-.#1)  due  to  an  overflow  In  the  format  print 
condition  of  the  computer  printout. 

The  correlation  matrix  is  a convenient  way  to  display  all  tabulated 
associations,  (R)  values,  whenever  two  or  more  seta  of  data  variables  are  being 
evaluated . 

A. 5 Correlation  Coefficients  of  Transformed  Data 

The  correlation  coefficient  (or  matrix  array)  is  a numeric  procedure 
whereby  the  potential  relationship  between  two  seta  of  data  can  be  examined. 
Nonlinear  as  well  aa  linear  relatiofnships  can  be  tested  for  their 
relative  degree  of  association  If  the  data  are  first  transformed.  Data 
Interdependence  which  can  be  linearised  on  log-log  or  semi-log  plots  can  be 
Inspected  analytically  by  taking  the  appropriate  logarithms  of  the  data  be- 
fore computing  the  correlation  coefficients.  Systems  governed  by  the  laws  of 
physics  can  often  be  explained  by  one  of  these  functional  relations.  Testing 
i for  one  of  these  relations  at  the  onset  of  examining  the  data  can  provide  gutek 

review  of  the  most  obvious  data  dependencies. 

A. 6 Selecting  an.  "Optimum  Set"  of  Explaining  Parameters 

The  "best"  regression  equation  Is  the  one  which  cen  sufficiently  describe 
the  system  under  study  with  the  amallest  nui^er  of  variable  parameters.  There 
Is  no  unique  statistical  procedure  that  will  provide  the  "beat"  descriptive 
equation  In  any  given  Instance.  Personal  judgment  Is  always  required.  Scien- 
tific Insight  and  engineering  principles  eheuld  never  be  substituted  by  some 
automatic  statistical  screening  method  for  obtaining  the  "best"  model  of 
representation  of  a given  set  of  data.  Some  analytical  procedures  are  more 
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uMful  than  othars  in  helping  find  an  "optlmun  set"  of  paraoeters  that 
can  sufficiently  describe  the  system  or  data  being  studied.  Two  types 
of  teats  are  cosnon;  parameter  acceptablliky  tests  and  optimisation  criteria 
for  gathering  the  "best"  set  of  parameters. 


Single  or  multiple  parameter  acceptability  can  be  established  for  a 
linear  set  of  data  through  the  use  of  the  F-test 

p ^ Mean  Square  of  the  Regression  Components 
Mean  Square  of  the  Error  Components 


Where  a regression  and  error  component  are  defined  &s  follows: 
Regression  components  * - Y) 

A 

Error  components  ■ (Y^  - 


Each  is  shown  graphically  for  a single  perameter  fit  in  Figure  k-'S.  This 
ratio  cam  be  calculated  sequentially  as  parameters  are  being  added  to  the 
model  or  collectively  to  all  parameters  in  the  model. 

The  (F)  number  for  parameter  acceptability  can  be  established  by 
using  any  of  a number  ot  prepared  tables  which  show  the  level  of  significance 
for  any  number  of  data  points  and  paraamters  in  the  model.  In  general,  the 
higher  the  F**nuid>er  the  better  chance  there  will  be  of  the  aukiel  successfully 
describing  some  additional  set  of  data  similar  to  that  from  which  the  model 
was  developed. 


Other  analytical  means  are  ofton  used  to  establish  an  "<vtlmum  set"  of 
parameters  from  those  available.  The  methods  comaMMily  used  Involve  the 
detenainatloo  of  all  errors  associated  with  the  data  and  the  ej^lainlng 
model.  The  errors  being  defined  as  the  differences  between  the  actual  data 
and  the  corresponding  predicted  ones. 


Three  calculating  procedures  or  criteria  are  helpful  in  determining 
whether  one  coaOtination  of  parametera  is  better  than  another.  Each  is 

82 


I 


calculated  aa  followa: 


2 Sub  of  Squaraa  of  the  Brrora  Y,^^l  ~ ^i^ 

^ Total  Sum  of  Squaraa  ^^1  “ 

A 2 

Sum  of  Squaraa  of  tha  Brrora  “ 'j) 

ISE.  - — " ■ ■ ■"  

k (N-k-1) 


r (Sum  of  Squaraa  of  tha  Brrora) 

“ Mean  aum  of  Squares  of 

Brrora  for  all  Poaalble 
Variables  in  Regression 

where  K la  the  number  of  data  points  and  k la  the  number  of  parameters 
entered  Into  the  regression  model.  Note  that  when  no  parameters  are  In 
regression  the  sum  of  the  squares  of  the  errors  Is  taken  as  the  total  sum 
of  the  squares. 


The  following  comments  are  pertinent  to  the  use  of  the  three  criteria: 


Criterion 


Ccaenents 

Does  not  account  for  the 
number  of  variables  In 
model. 


Value  Approached  as  Humber  of 
Variables  In  Model  Increases 

Rapidly  approaches  1 


MSE^ 


Change  in  often  very 
small  as  nuwer  of 


variables  in  model 
Increases 

May  Increase  slightly  upon 
adding  model  parameters 

Changes  in  MSB.  often  very 
wall  as  number  of 
variables  Increase 

Optimum  test  for  finding 
models  that  give  the 
smallest  squarsd  error 

Minimum  C|^  occurs  for  models 
containing  either  bias  or 
random  errors. 


Approaches  a minimum  (;the 
variance  of  the  errors) 


Approaches  k ^ 1 for 
"bMt”  set  of  variables* 
Ck  ■>  mmiber  of  perMetera 
In  model.) 


f -A 


Th«  conputttd  values  of  shown  graphically  In  Figure 

A-4.  In  each  caae  thoae  coadi>lnatlons  of  variables  providing  a value  of 
2 

RJ.  MSE^.  or  C|^  that  lies  nearest  the  solid  line  (the  left  boundary  of  the 
shaded  regions  shown)  will  likely  be  ''.he  optimum  comblnetlon  of  variables  for 
modelling.  The  final  selection,  of  uurse,  being  made  from  a scientific 
knowledge  of  the  system  being  studied. 

A. 7 A Large  Nuid>er  of  Independent  Variables 

The  "optlanjm"  selection  process  cannot  always  be  used  when  a large  number 
of  Independent  parameters  exist.  Since  each  potential  Independent  parameter 
can  be  Included  or  excluded  from  the  model  there  are  (2  ) possible  linear 

regression  equations  that  exist  with  k possible  parameters.  In  addition,  a 
few  perasMtera  nay  become  a large  collection  of  different  parameters  If  cross 
products,  logarithms,  or  other  typical  ways  of  modifying  the  data  are  to  be 
examined.  Even  with  the  availability  of  large  computers  running  all  possible 
regression  models  for  review  becoiBea  unwieldy.  In  the  case  of  a large  number 
of  available  parameters  other  procedures  must  be  employed. 

Among  tliose  used  ere  the: 

e)  Forward  selection  technique, 

b)  Backward  elimination  procedure,  and 

c)  Stepwise  procedure. 

The  objective  of  each  procedure  Is  to  give  the  best  description  of  the  response 
date  (dependent  variable)  with  the  fewest  number  of  convenient  peraeeters  In 
the  fuactlomel  expression.  Unfortunately  the  above  procedures  do  not  always 
lead  to  the  sene  solution  for  a given  set  of  data.  In  those  instesces  where 
the  date  variables  are  not  highly  interrelated  the  above  proceduree  will 
often  achieve  the  same  eolution.  In  cases  of  hl^^y  interrelated  parameters 
a unique  description  of  the  date  la  not  alwnya  avallahle.  In  these  Inwcances 
the  ’’boat"  eolution  will  require  a value  Judgment  or  minimum  errors  criteria 
toat  to  make  the  flmal  oelactlom. 
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The  forward  selection  technique  Ineerts  varlsbles  Into  the  regresel^i 
Mdel  vintll  a satisfactory  equation  is  achieved.  The  order  of  insertion  is 
usually  based  on  a calculation  of  the  fractional  correlation  coefficient  be> 
tween  the  existing  eodel  and  the  parameters  not  yet  entered.  The  parameter 
with  the  nlghest  coefficient  is  entered  next. 

The  backward  elimination  method  begins  with  all  variables  in  regression. 
Parameters  are  then  removed  aequentlally.  The  procedure  stops  when  the 
explaining  equation  is  sufficient  to  describe  the  data.  Partial  correlation 
coefficients  and  F-values  are  used  to  establish  which » if  any,  of  the 
parameters  should  be  taken  out  of  the  regression  equation  at  each  step. 

The  stepwise  p.*^jcedure  is  an  improved  version  of  the  forward  selection 
process.  The  Improvement  Involves  a rer-examlnatlon  of  each  parameter  in 
regression  after  one  is  entered.  A variable  previously  entered  may  be  found 
to  be  unnecessary  after  other  parameters  have  been  brought  into  the  BM>del. 

The  partial  F-test  is  applied  to  each  variable  after  each  parametar  is  entered 
into  regression.  If  any  variable  is  shown  to  be  nonsignificant  at  any  stage* 
it  la  immediately  removed  from  the  model.  The  entry-removal  procedure  con- 
tinues until  the  model  is  found  to  be  acceptable. 

A further  improvement  of  this  standard  technique  can  be  made  by  running 
the  test  (see  previous  discussion)  on  those  parameters  that  are  in  the 
model  when  the  stepwise  procedure  is  ended.  In  doing  so*  a fewer  nuid>er  of 
parameters  may  be  found  which  yield  an  equivalent  solution,.  In  any  case,  tl^s 
search  method  does  not  require  that  all  regresoions  be  analysed. 

A.  8 Regression  Tsblss  Bxplalned 

Every  useful  computerlied  regressioa  program  has  a talMlarlsed  output. 
Hoot  programs  provide  a similar  set  of  output  results  convenient  for  many  dif- 
ferent applications.  Figure  A-5  is  a sai^la  analytical  output  sheet  with 
annetated  explanation  of  tbs  pertinent  elements.  Horn  of  these  tables  will  be 
presented  below  wlttwxit  annotation. 
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A, 9 BxaapAes  of  Data  Analysis  Raduction 

As  shown  prsvlously  in  Flgura  I*  a knovlsdge  of  the  motor  temperature 
la  very  significant  In  estimating  how  long  It  %rill  last.  In  fact,  the  In* 
verse  of  the  absolute  motor  temperature  has  a correlation  coefficient  of 
R*.86  (See  Figure  A-2  column  16).  As  shown  In  Figure  A-2  this  Is  the  bast 
transfonsatlon  of  the  temperature  data  that  was  explored.  All  other  tempera- 
ture transformations  showed  a lower  degree  of  correlation  with  motor  life. 

The  practical  aspects  of  getting  available  temperature  data  tends  to 
make  the  information  provided  by  the  ambient  operating  temperature  more 
useful  than  the  others.  When  one  wants  to  predict  motor  lives  he  does  not 
always  have  the  motor  bearing  or  winding  temperature  at  hand.  The  ambient 
motor  operating  temperature,  however,  is  usually  available.  In  order  to  get 
a feel  for  the  relative  importance  of  ambient  and  bearing  temperature,  a 
single  variable  "fit"  to  the  life  data  was  performed  for  each. 

Tabular  outputs  of  the  regression  Information  for  the  bearing  related 
temperatures  and  the  ambient  operating  temperatures  are  shown  In  Figures 
A-8  and  A-7.  The  "best  fit  equations"  from  those  figures  are: 

Ambient  iMp.  Log  (G*  Life)  • 

Bearing  Temp.  Log  (C.  Life)  • 


2557 


W273 


- 2.534 


2342 


^273 


- 1.908 


Each  of  these  variables  are  shown  to  be  significant  In  estimating 
motor  life.  Both  have  F-values  In  emess  of  one  hundred. 

Bearing  tesqaratures  show  a better  "fit"  to  the  data  than  the  ambient 
(l.s.  f--value  of  159  versus  124).  in  addition,  once  the  bearing  temperatures 
arm  entarsd  Into  ragrssslon,  there  Is  no  Justification  for  using  ths  asdtlsnt 
or  winding  tsmperature  information.  This  conclusion  Is  drawn  from  the 
small  F<-values  related  to  the  sid>lent  and  winding  temperatures  shown  In 
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Figure  A-6  F-value  • .407  and  F-value 


.115  after  T.  had  been 
brg 

lf»  however,  the  ambient  tempera- 


accounted  for  in  the  regression  equation 
ture  has  been  entered  Into  regression,  there  Is  still  significant  ambient  or 
winding  temperature  Information  left  In  the  data  as  shown  In  Figure  A-7 
(T^d  F-value  ■ 10.0  and  F-value  ■ 18.7).  Even  though  the  bearing 
temperature  Is  more  significantly  correlated  with  the  life  data  as  a practical 
engineering  matter  It  was  decided  that  the  ambient  temperature  should  be 
used  in  estimating  motor  lives.  Since  ambient  motor  operating  temperatures 
are  more  readily  available  and  are  highly  significant  as  shown*  was 

chosen  as  one  of  the "primary  predictors"  of  motor  life. 


Other  parameters  chosen  as  life  predictors  were  reviewed  extensively 
In  the  main  body  of  this  report.  As  explained  in  the  previous  discussion, 
the  selection  of  an  optimum  set  of  "predictor"  variables  makes  use  of 
several  regression  procedures.  The  single  Interrelationship  between  the 
logarithm  of  motor  life  and  listed  variables  can  be  seen  In  Figures  A-8  and 
A-9.  Column  sixteen,  the  last  column  of  each  matrix  shown,  reveals  typical 
single  variable  dependence  (R  values)  with  respect  to  motor  life.  Some 
comments  are  appropriate. 

The  first  three  coefficients  In  column  sixteen  of  Figure  A-8  again 
show  the  dependence  of  life  on  absolute  temperature.  Temperature  shows  the 
strongest  correlation  of  any  of  the  variables  listed. 

The  fourth  eluent  of  column  sixteen,  however.  Is  misleading  as  shown. 
There  Is  no  apparent  dependence  (R-.OO)  of  life  on  the  D x N (bearing  diameter 
times  rotational  speed  of  the  shaft) . The  speed  dependence  does  not  show  up 
until  after  the  data  has  been  corrected  for  the  stronger  effects  such  as 
temperature  and  type  of  grease  used  In  the  motor  bearings. 

The  fact  that  motor  life  Is  dependent  on  D x N Is  shown  quite  well  In 
Figures  A-6  and  A-7.  Although  D x H Is  listed  as  "not  entered"  into  regres- 
sion, the  partial  F~values  are  shown.  These  partial  F-valuas  represent  the 
F-value  that  would  be  obtained  If  this  parameter  were  to  be  included  next 
In  the  explaining  equations.  The  F-values  of  3.8  (Figure  A-6)  and  5.4 
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(Figure  A-7)  represent  a highly  significant  dependence  between  D x N and 
motor  life.  Other  variables  found  to  be  significant  after  accounting  for 
temperature  were  the  quality  code  (listed  as  Itet;  7 in  Figures  A-7  and 
A-9)  and  grease. 


High  Negative  Correlation 


••  Y 
• 

a 

a. 


Some  Correlation 


Yl  • * • • 

v;> 


R&O 

fC>r  near  Zero) 


Figure  A-l 

Schematic  Plots  Depicting  Various  Aaounts 
of  Scatter,  as  Hsll  as  Correlation  Coefficient 
(R)  beteeen  Variable  X end  Response  Y 
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TEMPERATURE  CORRELATIONS  WITH  LOG (CHARACTER! STIC  LIFE) 


1 = AMBIENT 
5 = BEARING 

3 = WINDING 

4 1. /AMBIENT 

5 = 1 /BEARING 

6 = 1 /WINDING 

7 « 1/ (ABSOLUTE  AMBIENT) 

8 1/(.AB'3DIUTE  BEARING) 

9 * 1/ (ABSOLUTE  WINDING) 

10  - LOG (AMBIENT) 

IJ  = Lots  (BEAR  INC) 

15  = LOG(WINDING) 

13  - LOG(LOG(AMBIENT) ) 

14  = L{X(LO('(^EARING) ) 

15  = LOG ( LOG ( WINDING ) ) 

16  « I OG(CHARACT-ERiSTIC  LIFE) 

Figure  A- 2 

CORRELATION  MATRIX  FOR  OVER  1000  FHP  MOTOR  FAILURES 
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Si*  ■ 


N <■  Number  of  Data  Points 
1 B 1 th  Data  Point 


Figure  A'S 

Schamatlc  Data  Plot  Showing  Obaarvad  Data  Joints  (X^,  Y^), 
Bast  Fit  Pninta  (X^,  Y^),  Avaraga  x Y -Jy^/N 

Error  Valuaa  Naaldunls  (Y^-Y). 
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J<>  -<1  X> 


Ragr«88lon  Table  De8crlptlon 
93 


XXXXXXXXXXXXXXxXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
XXXXXXXXXXXXXXX  REGRESSION  TABLE  FOR  STEP  1 XXXXXX  T6RADC  XX 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
DESCRIPTION  TOTA_r  REGRESSION  RESIDUALS 


OBSERVATIONS 
CORRLATE  SS 
ORIC  UNITS 
CORR  MEAN  SOS 
GRIG  UNITS 
F TEST  NUMBER 


13S 


S835e0€.94 


1 

. 5435E3Ct00EBl 

25. 38699767946 
.S43S33800S51 

29.38699767946 
159.5S3D3230S1 


134 

.456476199749 
24.59656301454 
3. 40653880E--03 
. 183556440407 


XXXXXXXXXXXXXXX aXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


VARIABLES  ENTERED 
T BRG  ABS-1 


NOTE: - 


F-VALUES 
59. 5530923051 

VARIABLES  OUT  OF  REGRESSIaN^^\  F-VALUES 
" AMB  ^S-1  (NOT  ENTERED)  ^ 

WND  ABS-1  (NOT  ENTERED) 


1 AMB  MS-1 
T MND  ABS-1 
EARLY  »1 
LATE  -1 
U50  »1 
B35  =1 
033  =1 
0 X N 


(NOT 

(NOT 

(NOT 

(NOT 

(NOT 

(NOT 

(NOT 

(NOT 


4075622283483 

1153236849238 


ENTERED)  38.6195427391 

ENTERED)  38.61954272783 

ENTERED ) 42. 78808757997 

ENTERED)  .9609179303038 

ENTERED)  41 . 1 387084556 

ENTERED)  3.82090228193 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
PERCENT  VARIATION  EXPLAIICD  54 

STD.  DEV.  or  RESIDUALS  .42843487301 

XXXXXXXXXXXXXXXXXXX  EGUATION  OF  REGRESSION  XXXXXXXXXXXXXXXX 


CUT  OFF  VALUES 
3.29 

R2  UI/RESPQNSE 

1.39454206E-03 

3.95465495E-04 

. 1027208313472 

.102720831324 

.1111095972538 

3.27435920E-03 

. 1078384091908 

l.r'476936E-02 


VARIABLE 


SLOPE 


STD  DEVIATION 


CONSTANTS 


r BRG  A8S-1  2342.767986407  185.47127948  5.550021179435 


I 


I 


RESPONSE  VARIABLE 
LOG  C.LIFE 


TOTAL  CONSTANT 
-1.908089974872 


Sr 


Flgtire  A-6 


Regression  Table  for  Bearing  Temperature 
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xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
XXXXXXXXXXXXXXX  REGRESSIJN  TABLE  FOR  STEP  1 XXXX' X T6RADC  XX 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
DESCRIPTION  TOTALS  REGRESSION  RESIDUALS 


OBSERVATIONS 
CORRLATE  SS 
ORIC  UNITS 
COBR  MEAN  SOS 
ORIC  UNITS 
F TEST  NUMBER 


53. 88356069A 


.A80712SBaa08 
SS.9024BS43A71 
. 48071 SS8SS08 
25.90248943471 
124.0457719465 


134 

.519287717792 
27.98107125929 
3.87528147E-03 
. 20881 39e4€2ie 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


VARIABLES  ENTERED 
T AWB  ABS-1  note 

VARIABLES  OUT  OF  REGRESSI 
T WND  ABS-1  <NOT  ENTERED) 

T BRG  ABS-1  <NOT  ENTERED) 

EARLY  =1  <NOT  ENTERED) 

LATE  =1  INOT  ENTERED) 

U50  =1  (NOT  ENTERED) 

B35  »1  (NOT  El'nERED) 

D33  =1  (NOT  ENTERED) 

D X N (NOT  ENTERED) 


NOTE: »124 

ESSIof^s^F-V 

TERED) 

TERED)  >klB. 


F-VALUES 

►124.0457719466 

F-VALUES 

MO. 04881945347 

118.76455763472 

39.0887788064 

39.08877878905 

51.99320038274 

.2886183445159 

54.23319985223 


CUT  OFF  VALUES 
3.29 

R2  W/RESPONSE 

3.64786549E-02 

6.42060601E-02 

.1179526223525 

. 1 17952622312 

.1459482311328 

1.12444680E-03 

.1504147480364 

2.03118007E-02 


D X N (NOT  ENTERED)  5.414027820727  2.03118007E-02 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
PERCENT  VARIATION  EXPLAINED  48 

STD.  DEV.  OF  RESIDUALS  .45696166647 

XXXXXXXXXXXXXXXXXXX  EQUATION  (3F  REGRESSION  XXXXXXXXXXXXXXXX 
VARIABLE  SLOPE  STD  DEVIATION  CONSTANTS 


T AMB  A8S-1  2357.544811901 


211.67471 


6. 175985477438 


RESPONSE  VARIABLE 
LOG  C.LIFE 


TOTAL  CONSTANT 
-2. 534054272875 


Flg.«re  A-7 


Rcgiresilon  T«bl«  for  Aoblont  Teaptratur* 
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99 

.9€> 

.97 

.09 

. 08 

. 07 

. Ht> 

.08 

. (8 

.01 

. 13 

. 37- 

.04 

.75 

. 85 

.73 

9C' 

.9*9 

. 98 

.13 

.09 

. 09 

.47 

.03 

. 60- 

•.04 

. 16 

.40- 

. 18 

.73 

.83 

.77 

97 

. 98 

. i 1 

. 0'9 

.08 

. 45 

.Or2 

.64- 

.04 

. 15 

.40- 

.08 

. 74 

.8<3 

.75 

09 

.13 

.11- 

. 3 5- 

. 18- 

.03- 

.01 

.03- 

.05 

.31 

.09- 

. 13 

.0C-. 

.07 

.OO 

08 

.09 

.09  - 

.35 

.99 

.89 

.88- 

.06 

. 13 

.33- 

.85- 

. 04  - 

.0*9 

. 19 

.88 

.84 

07 

.09 

.0S■- 

. 18 

.89 

. 99 

.8t.- 

.05 

.11 

.35- 

.11- 

.08- 

.03 

. 19 

.88 

. 17 

AS 

.47 

.A5- 

.03 

.88 

.et. 

.00 

.80 

. 10- 

.08 

.05- 

.48 

.57 

.46 

.64 

08 

.03 

.08- 

.01- 

.08- 

. 05  - 

.00- 

.##- 

.00 

.04- 

.06 

.00- 

.00 

.06 

.01 

.08 

88 

.80 

.8A 

.03 

.13 

. 1 1 

.80- 

.OO- 

-.## 

.86 

. 19 

.85 

. 13 

.77 

. 17 

.39 

1 

O 

.04- 

.04- 

. 05 

^ ' 

.35 

. 10- 

.04 

.86 

.99- 

. 15- 

.07 

.03 

. 10 

. 18 

.00 

13 

. 18 

. 15 

. 31  - 

.85- 

-.11- 

.08- 

.06 

. 19 

.15 

.99 

. 13 

.09 

.13- 

.08- 

.08 

37 

.40 

.40 

.09- 

.04 

.08 

.05 

.00 

. 85- 

-.07 

. 13- 

.m 

.09 

.89- 

•.  10 

. 18 

04- 

. 18- 

.08- 

.13- 

.09- 

.03- 

-.48- 

.00 

. 13 

.0? 

.09 

.09- 

. 10- 

-.88- 

.44 

7S 

.73 

.74 

.<X> 

. 19 

. 19 

.57 

.oc. 

,77 

. 10 

.13 

.89- 

•.  10- 

•.## 

.40 

. 63 

85 

.83 

.80 

.07 

. 88 

.88 

.46 

.01 

. 17 

. 18- 

.08- 

. 10- 

-.88 

.40- 

.44 

73 

.77 

.75 

.00 

.84 

. 17 

. 64 

.08 

.39 

.00- 

.08 

. 18- 

.44 

.63 

.44- 

ROW  AND  COI.UMN  VARIABLES  ASSIGNED  AS  FOLLOWS 


1 = AMEit  ABSOl-UTE) 

a ^ 1/(TF£NP  BRG  ABSOLLHE) 

3 l/<  TLhP  WNG  ABSOI_UTE> 

A - 0*N  (DIA.  .X  RPM) 

5 NI.IMBER  ON  TEST 
G-.  = NUMECR  FAILED 

7 = QLIALITY  CODE 

8 = slope 

9 =-  TEfjT  CORR  COEF 
10  = F-TEST  VALUE 
n WIRE  SIZE 

18  NUMBER  MOTOJl  POLES 
13  ^ GREASE  CODE 
1^  = TEST  DATE 
15  * SERIES 

1€.  - LOG  ( CHARACTERISTIC  LIFE  ) 

CORRELATION  MATRIX  FOR  OVER  lOOO  FH^  MOTOR  FAILURES 

Figure  A-8 
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. A1 

. TO 

. It 
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.05 

.83 

.07 

.49 

.48 

.51 

.95 

.65 

.76 

.93- 

. #*t 

. 3'^ 

. AO 

. ::<o 

. 1 7 

O 

ifi 

.06 

.8A 

.07 

.49 

.48 

.51 

.95 

.68 

.75 

.33 

. 33 

.ttt? 

. AO 

. 30 

. 1 7 

. SO 

.06 

.8A 

.07 

.49 

.48 

.51 

.95 

.67 

- 75 

. A1 

. AO 

. AO 

.11 

. 08  - 

■ . 1 s - 

.OA 

. A7- 

.05 

.75- 

.00 

.55 

.45- 

. 16 

.31 

. :30 

. 3C 

. 30 

.11- 

. A6 

. 35- 

. 19 

.33 

-55 

.35- 

.01- 

.03 

. 35 

.54 

.49 

. It. 

. 1 7 

. 17- 

. 08 

. A6- 

■.m 

.35- 

.03 

.35 

.46- 

.00 

. 16- 

. 55 

.35 

.54 

.11 

. A3 

. SO 

. SC  - 

■.  IS 

.35 

. 35 

. 30 

. 17 

. 3A 

. 15- 

.06 

.50- 

.54 

.55 

.59 

.51 

.OS 

.06 

.06- 

.OA- 

•.  19- 

■ . 09 

. 17- 

.♦»#- 

.04- 

. 16- 

.09 

.45- 

.09 

. 14 

. 18- 

. 17 

.83 

. PA 

.8A 

. A7 

. 33 

. 35 

. 3A  - 

.OA- 

.34 

.55 

.36 

.59 

.89 

.60 

.60 

. 07 

. 07 

.07- 

. OS 

. 5S 

. A6 

. 15- 

■.  16 

.34- 

-.m- 

.OO- 

.05- 

-.  16 

. 18 

. 15 

. 14 

. A9 

. A3 

. A9 

. 7E 

. 35- 

. oo- 

-.06- 

.09 

.S5- 

-.00- 

.## 

.05 

.45 

.47 

.54 

.36 

, A8 

. A8 

. A8- 

.00 

. 01 

. 16 

. SO 

.AS 

. 36- 

.05 

.05 

.99- 

.03 

.53 

.56 

.05 

.51 

.51 

.51 

. S5  - 

-.03- 

. 5S 

.5A- 

-.09 

.59- 

■.  16 

.45- 

-.03 

.99 

.51 

.OO 

.13 

.95 

.95 

.95 

. A5 

. 35 

. 35 

.S5 

. lA 

.89 

. 18 

.47 

.53 

.51- 

.68 

.63 

.6S 

. 68 

.67- 

16 

.5A 

.5A 

.59 

. 18 

.ED 

. 15 

.54 

.56 

.00 

.68 

.99 

.40 

.76 

.7S 

.75 

. 31 

. A9 

. 11 

.51- 

•.  17 

.60 

. 14 

.36 

.05 

.13 

.63 

.40 

.99 

RLWCOLIJMI'J  VAR  ISLES  AS  FC1LL.0».>;S 

1 = TEMP  AMB  ABSt.lLUTE 
E ^ TEMP  DRG  ASSOI-UTE 
3 ==  TEMP  WNG  ABSC^LUTE 
A = D*N  <DIA.  X PPM) 

5 NIWBER  ON  lEST 
E NUMBER  FAILED 

7 ==  QUAL-ITY  CODE 

8 « I'JEIBUL  SLOPE 

9 TEST  CORR  COEF 
10  F-TEST  VALUE 
n = WIRE  SIZE 

15  ■-=  NUMBER  MOTOR  POLES 
13  " GREASE  CODE 
lA  TEST  DATE 

15  = SERIES 

16  = LOG  < CMARACTERISTIC  LIFE  ) 

CTOFtF^  fv|/s»'rF«XX  UlirMOIlMCS  I IL-UIReS 

Figur*  A-9 
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APPENDIX  B 


RECOMMENDED  REPLACEMENT 
OF 

RELIABILITY  MODEL  FOR  ELECTRONIC 
COOLING  DEVICES  AND  FHP  MOTORS 

FOR 

MIL-HDBK‘217B 

B.I  Motor  Reliability  Prediction  Model 

This  section  describes  the  method  to  be  used  for  calculating  failure 
rates  of  motors  with  power  ratings  below  one  horsepower.  The  types  of  motors 
that  the  predictive  techniques  were  developed  from  include  polyphase,  capacitor 
start  and  run  and  shaded  pole.  Its  application  may  be  extended  to  other  types 
of  fractional  horsepower  motors  utilising  rolling  element  grease  packed  bear- 
ings. The  model  is  dictated  by  two  failure  inodes,  i.e.;  bearing  failures 
and  winding  failures.  Application  of  the  model  to  D.C.  brush  motors  assumes 
that  brushes  are  inspected  and  replaced  and  are  not  a failure  mode.  Motors 
included  using  these  models  cover  applications  of  electronic  cooling  using 
fans  and  blowers  as  well  as  other  motor  applications.  Reference  (1*)  con- 
tains a more  comprehensive  treatment  of  motor  life  prediction  for  motor  appli- 
cations where  continuous  operation  at  extremes  of  temperature,  speed  or  load 
are  anticipated.  The  reference  should  be  reviewed  when  bearing  loads  exceed 
10  percent  of  rated  load,  spAads  exceed  24,000  rpm  or  motor  loadi.  include 
ontor  speed  slip  of  ^yreater  than  25  percent. 

Failure  rates  experienced  by  motors  are  not  constant  but  Increase  with 
time.  The  failure  rates  described  In  the  section  are  average  rates  obtained 
by  dividing  the  cumulative  hazard  rate  by  the  time  period  (t)  of  interest. 
Failure  rate  is  most  significantly  influenced  by  operating  temperature.  Care 
^Aihould  be  exercised  in  using  the  best  estimate  of  temperature  for  determination 
of  failure  rata. 


*5ee  Bibliography  Nuodter 
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The  failure  rate  nuxlcl  Is: 


t 


'sub  __1 
3 a 
B 


X 10^  (fallures/lO^hours) 


(1) 


where: 


the  failure  rate 


( fal lures/ 10 


6 


hours) 


■ tine  period  for  which  average  failure  rate  Is 
calculated  (hours) 


t^  * Estimated  average  failure  rate  at  MTBF  as 
tabulated  in  Table  A calculated  from 
Equation  (4).  Tables  are  used  for  constant 
temperature  operation  and  equations  (2), 

(3),  and  (4)  used  for  cyclic  temparsture 
operation. 

t2  * If  failure  rate  using  t.  is  excessive  relative 
to  equipment  requirements,  user  may  select  a 
smaller  value  of  t;  l.e.:  t..  All  motors  how- 

ever, must  be  replaced  at  t.  to  make  this  failure 
rate  valid.  User  selects  tnls  value,  it  Is  not 
calculated. 

Og  ■ Bearing  Welbull  Characteristic  Life  as  determined 
from  Table  A or  Equation  (2). 

o ■ Winding  Welbull  Characteristic  Life  as  determined 
* from  Table  A or  from  Equation  (3). 

B.2  Predicting  Failure  Rate  Without  Utilising  Tables 

Table  A tabulates  the  Welbull  characteristic  bearing  life  and  winding 

life  ( Ob  ■nd  a , respectively).  If  tables  are  not  used,  the  following  pro- 
B w 

cedure  may  be  followed: 


1.  Determine  operating  ambient  temperature  in  degrees  centigrade. 

Calculate  _ and  a from 

B V 


o 

B 


[2357 


10  > 


! 14-273 


2.5341 


+ ^-4500 

10  Tf273 


■h  20 


4-300 


(2) 
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MOTOR  R/Ml_tJRE  rate:  TABLE 

TABLE  A 


AMBIENT 

TEI*f»ERATURE 

<DEG  C) 

BEARING 

LIFE 

ALPI-tA  8 

WIMDING 
LIFE 
ALPHA  W 
(HRS) 

MOTOR  FAILURE  RATE 
IVRJMBER 

PER 

(MILLION  HOURS) 

140 

1489 

7533 

574 

138 

1587 

8031 

538 

136. 

1693 

8567 

504 

134 

1807 

9144 

472 

132 

1950 

9766 

442 

130 

2063 

10438 

414 

128 

2207 

11163 

387 

126. 

2361 

11947 

361 

124 

2529 

12794 

338 

122 

2710 

13711 

315 

120 

2907 

14704 

294 

118 

3119 

15780 

274 

116 

3350 

16948 

255 

114 

3601 

18215 

237 

112 

3873 

19591 

220 

no 

4169 

210tE(8 

205 

108 

4490 

22717 

190 

106 

4841 

24490 

176 

104 

5223 

26423 

163 

102 

5640 

28532 

151 

100 

6095 

30834 

140 

38 

6592 

33350 

129 

96 

7136 

36102 

119 

94 

7731 

39114 

no 

92 

8384 

42416 

101 

90 

9099 

46037 

93 

88 

9884 

50013 

86 

86 

10747 

54382 

79 

84 

11695 

59188 

73 

82 

12740 

64481 

67 

80 

13890 

70315 

61 

78 

1515P 

76753 

56 

76 

16555 

83865 

51 

74 

181v>4 

91729 

47 

72 

19812 

100434 

43 

70 

21700 

110082 

39 

68 

23790 

120787 

35 

66 

26101 

132679 

32 

64 

28658 

145903 

29 

62 

31484 

160627 

27 

60 

34605 

177042 

24 

<DEG  C) 


(MILLION  HOURS) 


(HRS^ 


(HRS) 


58 

38044 

195364 

88 

56 

41881 

815840 

80 

54 

45348 

838754 

18 

58 

50487 

864488 

16 

50 

55834 

833833 

15 

48 

60J17 

385596 

14 

46 

65571 

368005 

18 

44 

70887 

403084 

11 

48 

75889 

449303 

11 

40 

80884 

501598 

10 

38 

83575 

560761 

10 

36 

85338 

687815 

9 

34 

85856 

703981 

9 

38 

88875 

790439 

9 

30 

78855 

888951 

10 

88 

71716 

1001301 

11 

86 

63838 

11896.48 

18 

84 

55318 

1876517 

14 

88 

46818 

1444875 

16 

80 

38887 

1638806 

SO 

18 

31677 

1860614 

84 

16 

85508 

8116943 

30 

14 

80316 

8418983 

38 

18 

16050 

8755341 

48 

10 

18538 

3158859 

68 

8 

9840 

3613869 

79 

6 

7657 

4149814 

108 

4 

5944 

4775570 

131 

8 

4609 

5506984 

169 

0 

3574 

6363558 

819 

-8 

8776 

7369148 

888 

-4 

8165 

8558877 

361 

-6 

1698 

9947538 

461 

-8 

1343 

11596838 

583 

-10 

1075 

13551188 

788 

-18 

873 

15878719 

897 

-14 

788 

18637433 

1084 

-16 

609 

81938464 

1886 

-18 

585 

85890310 

1491 

-80 

463 

30634188 

1689 

-88 

417 

36344468 

1875 

-84 

384 

43837786 

8036 

-86 

360 

51583418 

8174 

-88 

348 

61717476 

8885 

-30 

330 

74060813 

8370 

-38 

381 

89148185 

8439 

-34 

314 

107687855 

8489 

-36 

3J0 

130360998 

8585 

-38 

30? 

158411698 

8549 

-40 

304 

193143380 

8570 

101 

(3) 


<*  . ^ . 1 83 

w Tf273 
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2.  Utilising  the  above  equations  fora  and  a , the  mean  time  between 

B W 

failure  (t)  may  be  estimated  from; 


t ■ o. 


N 34657  34657)^  + . 03704 f 


/. 34657  -\/(.3A657V 


+ . 03704  r 

a 

L wJ 


3.  Calculate  the  failure  rate  from  Equation  (1), 


B.2.1  Bearing  Characteristic  Life  Resulting  from  Thermal  Cycling 


(hj  4 hg  + h3  hj 

^ **1  . . **3  I ■ ■ - h 

“•i  “-a  % t 


idtere:  (See  Figure  9-1) 


• time  at  temperature 

" tlsM  to  cycle  from  tmaperature  to 

"•  time  at  temperature 

• time  at  tmeperature  Tg| 

bearing  life  at  from  Table  A 

• bearing  life  at  T2  from  Table  A 
" bearing  life  at  from  Table  A 


‘In  r ■ '-Jj  • 


. ■■ 


I . ^*^3 

2 2 


from  *xampl«  of  Figure  A-l, 


■■‘'4 — 2 — i*"^  \ ■ \ 


B.3  Winding  Characteristic  Life  Resulting  from  Thermal  Cycling 


(h.  + h-  + h,  + - 
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where:  (See  Figure  A-)) 


time  at  temperature 

time  to  cycle  from  temperature  to 

time  at  temperature 

time  at  tmaperature 

m 

winding  life  at  from  Table  A 
winding  life  at  T^  to  T^  from  Table  A 


winding  life  at  T from  Table  A 

ID 

Tj  + Tj 


from  example  of  Figure  A-1  T^^  - T^  - ^3  ^1  and  h^^ 


B,4  Examples  of  Kotor  Failure  Rate  Prediction 


B.4.1  Example  Wo.  1.  Given:  Fractional  Horsepower  Motor  operating  at  an 
agblent  tamperature  of  120°C.  Find  fallnre  rate  for  replacement  of  motors 


t ,•  :«•, , « 
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In  1000  hour*.  Also  find  fallurn  rate  at  eatlaiated  naan  tine  between  fall- 


a)  For  replacement  of  motors  In  1000  hours; 
Froei  Table  A at  120°C 


2907  hours 


* 14704  hours 


From  k^atlon  (1) 


t a J a 

B w 


jml  + -1-  X 10® 

2907’ 


" 108.7  failures/ 10  hours 


b)  For  failure  rate  at  estimated  MTBF 
From  Table  1 at  120°C 


■ 294  failures/ 10  hours 


en;  Fractional  Horsepower  Motor  operating  at  a 


thermal  duty  cycle  of: 


2 hours  at  100*’c  aadilent 

8 hours  at  20**C  ambient 

O.S  hours  to  change  temperature 
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Find  average  failure  rate  at  MTBF 

a)  Tabulata  operating  temperatures: 


Tj  ■ 100  C;  hj^  • 2 hours 


100  + 20  ,-0_  . , . 
■ 60  C;  ^2  ^ 


T«  ■ 20  C;  h_  ■ 8 hours 


b)  Detemine  bearing  and  winding  life  from  Table  A 

Tj  ■ lOO^C;  Og  “ 6095  hours;  ■ 30834  hours 


'^2  ■ 60  C;  “g  “ 34605  hours;  ci^  *•  177042  hours 


T_  ■ 20^C;  a_  ■ 38827  hours;  • 1638206  hours 


c)  Calculate  bearing  life  and  winding  life  from  Equations  (5) 
and  (6) 


“B 


2 + 1 + 8 

2 1 8 
6095  34605  38827 


19535  hours 


a 

w 


2 + 1+8 

2 1 8 
30834  177042  1638206 


145899  hours 


d)  Calculate  (t)  from  Equation  (4) 
t - 19535 


"V 34657  +"\/^201  + 03^  yHHo  ^ + 


/s^es?  - V*12  + .03^  mil. 


<1  J 

- 19535  V*  69325  + \-.000ll 

- 19535  ^.88504  - .0483] 

t *■  16346  hours 
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e)  Calculate  average  failure  rate  at  MTBF  from  Equation  <1) 
2 


\ " 


a 

V 


X 10 


(16346)  . 1 

(19535)3  '”899 


X 10‘ 


■ 42.7  failures/ 10°  hrs 


